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Abstract: Recently, an explicit, recursive formula for the all-loop integrand of pla- 
nar scattering amplitudes in M=A SYM has been described, generalizing the BCFW 
formula for tree amplitudes, and making manifest the Yangian symmetry of the the- 
ory. This has made it possible to easily study the structure of multi-loop amplitudes 
in the theory. In this paper we describe a remarkable fact revealed by these investi- 
gations: the integrand can be expressed in an amazingly simple and manifestly local 
form when represented in momentum-twistor space using a set of chiral integrals 
with unit leading singularities. As examples, we present very-concise expressions for 
all 2- and 3-loop MHV integrands, as well as all 2-loop NMHV integrands. We also 
describe a natural set of manifestly IR-fmite integrals that can be used to express 
IR-safe objects such as the ratio function. Along the way we give a pedagogical 
introduction to the foundations of the subject. The new local forms of the integrand 
are closely connected to leading singularities — matching only a small subset of all 
leading singularities remarkably suffices to determine the full integrand. These re- 
sults strongly suggest the existence of a theory for the integrand directly yielding 
these local expressions, allowing for a more direct understanding of the emergence 
of local spacetime physics. 
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1. Invitation to Local Loop Integrals and Integrands 

The integrand for scattering amplitudes in planar theories is a well-defined, rational 
function of external- and loop-momenta at all orders of perturbation theory [1]. 
Recently, an explicit recursion for the integrand of planar scattering amplitudes in 
M = 4 SYM was presented [1], generalizing the BCFW recursion for tree amplitudes 
[2,3]. The integrand is most naturally presented in momentum-twistor space. All the 
objects appearing in the recursion relation have simple interpretations in terms of 
canonical operations on Yangian-invariants derived from the Grassmannian integral 
[4], making the Yangian invariance of the theory (up to total derivatives) manifest at 
the level of the integrand. It has also been recently realized that the integrand has 
a beautiful dual interpretation as a natural supersymmetric Wilson loop, resolving 
a long-standing open problem [5,6]. This proposal has been checked to satisfy the 
all-loop recursion relation at the level of the integrand [5], providing a proof of the 
duality between scattering amplitudes and Wilson-loops [7]. 

The recursion relation gives a complete definition for the integrand, making no 
explicit reference to spacetime notions either in the usual or dual spacetimes. The 
words "spacetime" , "Lagrangian" , "path integral" and "gauge symmetry" make no 
appearance. A reflection of this fact is that, as familiar from the BCFW computation 
of tree amplitudes, individual terms in the integrand are riddled with non-local poles 
that cancel in the sum. But also familiar from BCFW at tree-level, the recursion 
relation is a very powerful calculational tool, and has allowed us to gather a huge 
amount of "data" about the properties of multi-loop amplitudes. 

In this paper we report on a remarkable property of the loop integrand revealed 
by examining this "data" , amplifying a theme already stressed in [1]. Loop integrands 
take an amazingly simple form when expressed in a manifestly local way. This is 
surprising, since the enormous complexity of Feynman diagrams is inexorably tied to 
locality, while by contrast, the great simplicity of BCFW recursion is inexorably tied 
to the presence of non-local poles. What we are finding is a new local form of the 
integrand — certainly not following from Feynman diagrams! — which is even simpler 
than the forms obtained from BCFW recursion. 

This great simplicity is apparent only when the integrand is written in momentum- 
twistor space, using a special set of objects that are almost completely chiral, and 
have unit leading singularities. For instance, all 2-loop MHV amplitudes are given 
as a sum over a single type of object, 
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This result was already presented (albeit in a slightly more clumsy form) in [1]. We 
will describe these objects in much more detail in the body of the paper; here, it 
suffices to say that these are simple double-pentagon integrals with a special tensor- 
numerator structure which is indicated by the wavy lines, and that the notation 
l i<j<---<k<i' in the summand should be understood as the sum of all cyclically- 
ordered sets of labels i,j,...,k for each % e {1, . . . , n}. 

All 2-loop NMHV amplitudes are also associated with similar integrands; indeed, 
the n-point NMHV scattering amplitude's integrand is simply given by, 




Here, [ijklm] denotes the familiar dual-superconformal invariant of five particles, 



\i ' k I ml = ^° |4 ((jklm)rji + (klmi)^ + (lmij)rj k + (mij k)rji + {ij kl)rj m ) 
^ (ijk I) {j k I m) (k I mi) (I mi j)(mi j k) 

(1.3) 

This result dramatically simplifies the way this result was presented in [1] for the fr- 
aud 7-particle 2-loop NMHV integrands. 

Finally, all 3-loop MHV amplitude integrands are given by a sum over the same 
types of objects, 
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These explicitly-local, manifestly cyclic results for all 2-loop NMHV and 3-loop 
MHV amplitudes are new, and stunningly-simple — even simpler than the form pro- 
duced by the loop-level recursion formula. 

As we will see, these extremely simple expressions are very closely related to the 
leading singularity structure of the theory. The reason for the dramatic simplicity 
of these results relative to the ones presented in [1] is that in [1], each integrand 
was straightforwardly expanded in terms of a fixed basis of chiral integrals with unit 
leading singularities, while here we are tailoring the objects that appear directly 
to the amplitude. The structures are motivated by matching a particularly simple 
set of leading singularities of the theory; this is made possible only by using chiral 
integrands with unit leading singularities, which is why these objects play such a 
crucial role in the story. What is remarkable is that matching only a small subset 
of leading singularities in this way suffices to determine the full result. Of course, 
we confirm this not by laboriously matching all leading singularities, but rather by 
directly checking the conjectured local forms against what we obtain from the all-loop 
recursion relation. 

We do not yet have a satisfactory understanding for the origin of this amazing 
simplicity. Certainly, these expressions differ from the BCFW form in that they 
are not term- by-term Yangian invariant. This suggests the existence of a deeper 
theory for the integrand that will directly produce these new local forms, allowing a 
more direct understanding of the emergence of local spacetime physics. We strongly 
suspect that it is this formulation that will also help explain the amazing simplicity 
[9] seen in the integrals yielding the physical amplitudes, and also form the point 
of contact with the remarkable integrable structures of M = 4 SYM — Y-systems 
and Yang- Yang equations — seen at strong coupling and also in some collinear limits 
[10-12]. 

In [9], a geometric picture for scattering amplitudes is advanced, building on a 
beautiful paper of Hodges [13], which may shed some light on the origin of these 
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new local expressions. Hodges interpreted NMHV tree amplitudes as the volume of 
certain polytopes in momentum-twistor space, and showed that a natural class of 
triangulations of this polytope correspond to different BCFW representations of the 
amplitude. In [9], it is shown that at an even simpler triangulation of the same poly- 
tope is possible, yielding a new, manifestly-local formula for NMHV tree-amplitudes. 
Also in [9], a completely analogous 'polytope' formulation is presented for all 1-loop 
MHV amplitudes. Again, one natural set of triangulations leads to the BCFW form 
of the integrand, while even simpler triangulations directly lead to a number of new, 
manifestly local forms for the integrand. While this polytope picture has not yet 
been generalized beyond these most elementary cases of NMHV tree and MHV 1- 
loop amplitudes, the extremely simple local forms for higher loop amplitudes we 
present in this paper strongly encourages the thought that an appropriate extension 
of this idea must be possible. 

We should stress that when we say our results for the integrand are "manifestly 
local", we mean that the poles involving the loop integration variables are local. Of 
course the integrand should be "ultralocal" , that is, the poles involving both the 
loop integration variables as well as the external momenta must be local. The MHV 
integrands we present trivially have this property, but for NMHV amplitudes, our 
expressions involve the standard i?-invariants which have spurious poles as function 
of the external particle momenta. Given the beautiful, local form of the NMHV tree 
amplitude obtained from the polytope picture [9], it is quite likely that there is an 
even nicer representation of loop amplitudes which are not only local but ultralocal. 
This fascinating possibility certainly merits further exploration, but is beyond the 
scope of the present paper. 

We close this invitation with an outline for the rest of the paper. We begin with 
a pedagogical introduction to some of the foundations of the subject in section 2 
starting with a review of momentum-twistors and some of the associated projec- 
tive geometry in CP 3 . We also discuss how planar loop integrals are written in 
momentum-twistor space; while our focus in the paper is on M = 4 SYM, we ex- 
pect that the momentum-twistor representation of loop amplitudes will be extremely 
useful for any planar theory. We discuss the way that momentum-twistors make in- 
tegral reduction trivial, and illustrate this by showing how the 1-loop integrand can 
be reduced to a sum over pentagon integrals. Finally we discuss leading singular- 
ities at 1-loop and beyond in momentum-twistor language. The standard exercise 
of determining quadruple-cuts in momentum space is mapped in momentum-twistor 
language to a simple, beautiful and classic problem in enumerative geometry first 
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posed by Schubert in the 1870's, and we discuss the solution of these "Schubert 
problems" in detail. 

In section 3 we introduce chiral integrals with unit leading singularities which 
play a central role in our story. We illustrate how they work starting with the 
simplest case of 1-loop MHV amplitudes. 

In section 4, we discuss another feature of chiral integrals with unit leading 
singularities — generic integrals of this form are manifestly infrared finite, and can 
be used to express finite objects related to scattering amplitudes, such as the ratio 
function [14]. 

In section 5, we construct a basis for all 1-loop integrands, whose building blocks 
are not the familiar boxes or even pentagons, but a natural set of chiral octagons 
with unit leading singularities. We also compute the finite 1-loop integrals explicitly, 
and use these results to give a simple formula for the NMHV ratio-function at 1-loop, 
for any number of particles. 

In section 6, we discuss multi-loop amplitudes. We describe our heuristic strategy 
for using leading singularities to tailor momentum-twistor integrals to the amplitude, 
and show how this works for the 1-loop MHV amplitude, reproducing one of the 
local forms first derived using the polytope picture of [9] . We also discuss the 1-loop 
NMHV amplitudes in the same way. We then extend these methods to two loops 
and beyond, and show how to "glue" the 1-loop expressions together to produce 
natural conjectures for all 2- and 3-loop MHV amplitudes, as well all 2-loop NMHV 
amplitudes. These conjectures are verified by comparing with the integrand derived 
from the all-loop recursion relation. 

A number of appendices discuss various technical points needed in the body of 
the paper, including a detailed discussion of the 2-loop NMHV and 3-loop MHV 
integrands. 

2. Foundations 

In theories with massless particles, a well-known and convenient way of trivializing 
the constraint p\ = for each particle is to introduce a pair of spinors X^ and \( a \ 
replacing i->- (p a )aa = Pa{ a n)aa = Aa^A^. Of course, this map is not invertible, as 
any rescaling {A, A} — > {tX, t~ 1 X} leaves p invariant. This reflects that these variables 
come with a new source of redundancy; in the case of particles with spin, this re- 
dundancy is quite welcomed as it allows the construction of functions that transform 
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with fixed projective weights as S-matrix elements under Lorentz transformations. 
This is all well-known under the name of the spinor-helicity formalism [15-19]. 

Amplitudes are supported on momenta that satisfy momentum conservation. 
Clearly, it would be convenient to find variables where this constraint, ^ a p a = 0, is 
trivial. In planar theories, where color ordering is available, there is a natural way to 
achieve this, by choosing instead to express the external momenta in terms of what 
are known as dual-space coordinates, writing p a = x a — x a -i, [20]. 

To see the role played by planarity, consider the standard decomposition of scat- 
tering amplitudes according to the overall color structure, keeping only the leading 
color part: 

A n = Tr(T ai T a2 . . . T an ) A n (l, 2, . . . , n) + permutations; (2.1) 

here, each partial amplitude A n {l, 2, . . . , n) can be expanded in perturbation theory, 
and we denote the L-loop contribution by A^~ loop . Partial amplitudes are computed 
by summing over Feynman diagrams with a given color-ordering structure. 

In this paper we only consider the planar sector of the theory, and therefore 
A^~ loop will always refer to the leading-color, partial amplitude in the planar limit. 

Restricted to a particular partial amplitude, say, A n (l, 2, . . . , n), each momenta 
can be expressed as the difference of two "spacetime" points. More precisely, we 
make the identification p a = x a — x a -i, with p\ — x\ — x n . It is clear that mo- 
menta obtained in this way automatically satisfy ^2 a p a — — and the redundancy 
introduced in this case is a translation x a — > x a + y by any fixed vector y. 

Now, the only poles that can occur in A n (l, 2, . . . , n) are of the form Ylm=aPmi 
i.e., only the sum over consecutive momenta can appear. In the dual variables these 
become Ylm=a+iP m = Xa ~ Xb - The same kind of simplifications happen in planar 
Feynman diagrams to all orders in perturbation theory as we will describe. 

Now we have the variables {A, A} which make the null condition trivial while ig- 
noring momentum conservation, while the dual-space variables do the opposite. It is 
perfectly natural to wonder if there exists any way to combine these two constructions 
which makes both the null-condition and momentum conservation trivial. It turns 
out that such a set of variables does exist: they are known as momentum-twistors 
and were introduced by Hodges in [13]. 

The standard twistor construction developed in the 1960's [21] starts by making 
a connection between points in an auxiliary space — twistor-space — and null rays 
in spacetime. Likewise, a complex line in twistor space is related to a point in 
spacetime. The key formula is called the incidence relation, according to which a 
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Figure 1: Defining the connections between momentum- twistors, dual-coordinates, and 
cyclically-ordered external four-momenta 




point x in spacetime corresponds to set of twistors Z = (A, //) which satisfy 

= x^aX-. (2.2) 

Twistors satisfying this relation form a projective line in CP 3 . Even though Z has 
the components of a point in C 4 , the incidence relation cannot distinguish Z from 
tZ, and therefore the space is projectivized. 

In order to specify a line in twistor space — and therefore a point in spacetime — 
all that is needed is a pair of twistors, say Za and Zb, that belong to the line. Given 
the twistors, the line or spacetime point is found by solving the four equations coming 
from imposing the incidence relation for Za and Zb with x. It is easy to check that 
the solution is, 

(Here, we have made use of the familiar Lorentz-invariant contraction of two spinors 
(Aa As) = e^/jA^A^). 

Hodges' construction starts with any set of n twistors {Z\, . . . , Z n }. Using the 
association x a 4-» (Z a , Z a+ i), n spacetime points are defined. Quite nicely, it is trivial 
that p 2 . = (x a — x a ^i) 2 = because the corresponding lines, or (CP 1 s), intersect. This 
is illustrated in Figure 1. 

Given the importance of this latter fact, it is worth giving it a slightly more 
detailed discussion than we have so far. If two lines in twistor-space intersect, i.e. 
share a twistor Z int , then the corresponding spacetime points, say x and y, associated 
with the lines are null-separated. To see this, take the difference of the incidence 
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relations for Z\ 




yaa\nti 



to get 



0; 



which means that the 2 x 2-matrix (x — y) has a non-vanishing null eigenvector, 
i.e. A^ t , and therefore the determinant of (x — y) vanishes. But the determinant is 
proportional to (x — y) 2 when x and y are taken as vectors; and therefore x and y 
are null separated. 

As useful background for the rest of the paper let us discuss the null-separation 
condition, which is a conformally invariant statement, in twistor space. Consider 
again two generic spacetime points x and y and choose two representatives of the 
lines associated to them in twistor space, say, (Za, Zb) and [Zc, Zd). Treating each 
twistor as a vector in C 4 there is a natural SX(4) (conformal) invariant that can 
be constructed. This is done by contracting all four twistors with the completely 
antisymmetric tensor €jjkl to produce 



Clearly, this conformally-invariant quantity must encode information about how 
x and y are causally related. The Lorentz invariant separation (x — y) 2 is not 
conformally-invariant because it is not a cross ratio. However, the way to relate 
the two quantities is simple 



This relation is consistent with our earlier finding that if the points x and y are 
null-separated, then the twistors Za, Zb, Zc and Zd, are coplanar as points in CP 3 . 
In other words, the two complex lines intersect. 

When twistors are used to produce a configuration of points in spacetime which 
are pairwise null separated and then used to build momenta, the corresponding 
twistor space is called momentum-twistor space [13]. 

This twistor construction is in fact slightly more involved when one is interested 
in real slices of spacetime. In our discussion so far, we have been assuming that 
momenta are complex and hence the dual spacetime is complexified. This is useful 
for e.g. defining the usual unitarity cuts of loop amplitudes. In this paper, the 
complex version suffices and we refer the interested reader to [13,22]. 



{ZaZbZcZd) — £ijklZ a Z b Zc Z 



(2.4) 




(2.5) 
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A related construction is called dual momentum twistor space. Here 'dual' refers 
to the usual geometric — 'Poincare' — dual of a space. In other words, the dual space is 
the space of planes in CP 3 . Points in the new space which is also a CP 3 are denoted 
by Wi. The construction maps points to planes and lines to lines. In Hodges' 
construction [13], there is a natural definition of dual points associated to the planes 
defined by consecutive lines of the polygon in momentum twistor space of Figure 1. 

The construction defines a dual polygon by introducing dual momentum twistors 
W a defined by 

(W a )i = tIJKLZ *- lZ * Z a+ 1 (2 6) 

This definition is made so that W a contains A a as two of its components. 
2.1 Loop Integrals 

The focus of this paper is loop integrands and integrals. Here too, it is well known 
that in planar theories, loop integrals are very naturally expressed in terms of dual 
spacetime coordinates. Consider a very simple 1-loop integral, known as a zero mass 
integral, 



w 



d L L 2 (L - pi) 2 (L - pi - P2) 2 (L -pi-p 2 - P3) 2 



where the external momentum at each of the four vertices is null (hence the name) 
and N = G°i+P2) 2 (P2+Ps) 2 is a convenient normalization factor. Momentum conser- 
vation gives Pi = —p\—pi —p-s, and introducing the dual-coordinates p a = x a — x a _i, 
it is easy to see that the unique choice of L that makes translation invariance (in 
x-space) manifest is L = x — x^. The integral becomes [20] 



'x x 



(x — Xi) 2 (x — X2) 2 {x — X^) 2 {x — X4) 2 ' 
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where N = (x\ — x^) 2 {x2 — X4) 2 . Imposing translation-invariance gives rise to the 
same integral in x-space regardless of the original definition of L in the loop diagram. 
In other words, a different propagator could have been chosen to be L and the form 
(2.8) would still be the same. This uniqueness plays a crucial role in the definition 
of the integrand of the theory. 

Integrating over all points x in spacetime is the same as integrating over all CP^s 
in CP 3 . As before, each line in twistor space can be represented by a pair of twistors 
x <H> (Za, Zb)- Clearly, any GL(2, C) transformation on the A, B "indices" leaves the 
line invariant. Therefore the integral over spacetime is the same as the integral over 
the pairs (Za,Zb) modulo GL(2). This is nothing but the Grassmannian G(2,4) 
which can be parameterized by a 2 x 4 matrix 



We can immediately write a measure which is GL(2)-invariant by integrating over 
all Za's and Z B 's together with a combination of 2 x 2 minors of the matrix (2.9) 
with total weight —4. It turns out that the precise measure that corresponds to a 
d A x integration is 



where (A^ A#) is the (1 2) minor of (2.9) — the determinant of the first two columns of 
the 2x4 matrix (2.9). In the twistor literature this is written as (Xa \b) = (ZaZb loo) 
where (Ioo) KL is the infinity twistor which is block diagonal with the only nonzero 
diagonal element equal to e a b- loo is called the infinity twistor because it corresponds 
to a choice of the point at infinity in spacetime and therefore a line in twistor space. 
Its presence therefore breaks conformal invariance. This is not surprising as the 
measure d 4 x 'knows about' the metric in spacetime. 

Since the integration over lines will appear in many different contexts in this 
paper we introduce a special notation for it. Let's define 



The reason we have not included the factor (A^ A^) 4 in the definition is that in this 
paper we mostly deal with Af = 4 SYM and in its integrand factors with infinity 
twistors cancel. 




(2.9) 




(2.10) 




(2.11) 



AB 



-11- 



Going back to the loop integral in x-space (2.8), one can introduce the four mo- 
mentum twistors in Hodges' construction {Z 1 , Z 2 , Z 3 , Z 4 } to describe the external 
particles. Using the relation between the Lorentz invariant separations and momen- 
tum twistor invariants in (2.5), the integral (2.8) becomes 



where (ijkl) stands for the determinant of the 4x4 matrix with columns given by 
four twistors Z^ Zj, Zj, Zk defined in (2.4). 

One of the remarkable facts about (2.12) is that all factors involving the infinity 
twistor have disappeared. This means that the integral is formally conformal invari- 
ant under the conformal group that acts on the dual spacetime. This is why it is 
said to be dual conformally invariant (DCI). 

Clearly, if we had started with a triangle integral then the factor (Z1IZ2) = 
(Ai A2) would not have canceled and would have remained with power one in the 
denominator as if it were a propagator. Indeed, this viewpoint trivializes the sur- 
prising connections made in the past between the explicit form of triangle and box 
integrals. In other words, one can think of a triangle integral as a box where one of 
the points is at infinity. 

Once again, a careful definition of the contour which should correspond to only 
points in a real slice of complexified spacetime is not needed in this paper. It suffices 
to say that on the physical contour, the integrals can have infrared divergences (IR). 
This is the reason why we said that the integral was 'formally' DCI. We postpone a 
more detailed discussion of IR divergences to section 4. 

The purpose of this section is to show how momentum twistors are the most 
natural set of variables to work with loop amplitudes in planar theories. In order to 
do this we will first show how many familiar results can be translated into momentum 
twistors. Not infrequently, momentum twistors will completely clarify physics points 
which have been misunderstood in the literature. 




(2.12) 



AB 
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Integral Reduction at One-Loop Level 

In a general theory, 1-loop integral reduction techniques allow scattering amplitudes 
to be expressed as linear combinations of a basic set of scalar integrals 1 . The integrals 
have the topology of bubbles, triangles or boxes. 

Let us start this section by translating each of the integrals in the standard basis 
into momentum twistor language. Their corresponding form in momentum twistor 
space is 



i + l 



'Box 



z + i- 



<i »+l>0' J'+l><* «:+!)<( i+l> 



AB 



(AB % i+l)(AB j j+l)(AB k k+l)(AB I 1+1} 



' Triangle 



'Bubble 




7+1 



i + l 



(i i+l){jj+l)(k k+1) 



AB 



(AB) (AB i i+l)(AB j j+l)(AB k k+1) 



(2.13) 



(i i+l)(j ]+l) 



AB 



(AB) 2 (AB i i+l)(AB j j+1)' 



Note that here we have translated the plain scalar integrals without any normaliza- 
tion factors. Once again, only boxes are dual conformal invariant except for an overall 
factor which only depends on the external data. This factor involving 2-brackets and 
hence the infinity twistor can always be removed by a proper normalization as done 
in the zero-mass example (2.12). Scalar boxes in momentum twistor space have also 
been recently studied in [22,24]. 

A well known fact about M = 4 SYM is that at 1-loop level, bubbles and triangles 
are absent and all one needs are scalar box integrals. However, as we will see, this 
point of view is not the most natural one and actually turns out to be misleading. 

In order to understand this point, one needs to review the reduction procedures 
used to reach this conclusion. Before doing that let us mention some useful facts 
about momentum twistors. 

lr rhis is true in theories with no rational terms or in general theories for what is known as the 
cut-constructible part of them. See [23] for more details. In AT — 4 SYM rational terms are absent. 
This is why we do not elaborate more on this point. 
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In loop integrals, combinations of momentum twistors of the form Z^Zg make 
an appearance in every expression (where the brackets mean that the indices are anti- 
symmetrized) , reflecting the fact that it is the line (AB) that is being integrated-over, 
and not the individual twistors Za and Z B . 

These two-index objects are a class of more general ones called bitwistors. A 
generic bitwistor is a rank-two antisymmetric tensor Y IJ . Given two bitwistors, Y 
and Y, the conformally-invariant inner-product is given by (Y Y) = €uklY ij Y kl . 
A bitwistor which can be written in terms of two twistors cis 2^ ^ ]} is called simple. 
It is easy to show that a bitwistor is simple if and only if Y 2 = with the product 
defined as above. 

The reason for discussing bitwistors is that they provide a very natural integral 
reduction procedure. The procedure can be applied to integrals at any loop order 
but in this section we concentrate on only 1-loop integrals. The procedure we are 
about to present is in part the momentum twistor analog of the one introduced by 
van-Neerven and Vermaseren in [25]. 

At one loop one starts with general Feynman integrals of the form 

where the tensor T is made out of polarization vectors, momenta of external particles 
and the spacetime metric. 

By Lorentz invariance, it is clear that one can decompose integrals of this type 
as linear combinations of momentum twistor tensor integrals of the form 

1 (ABY 1 )(ABY 2 )...(ABY m ) 

(ABI 00 )*-( n ~ m ') (AB 12} (AB 23} ...{ABn-1 n){ABnl) K ' ' 

AB 

where Y a are generic bitwistors. 

The reduction procedure relies on the fact that a generic bitwistor has six degrees 
of freedom and can therefore be expanded in a basis of any six independent bitwistors. 
To reduce the integrals in (2.15) simply choose any six of the bitwistors that appear 
in the denominator, say, Z\Z 2 , Z 2 Z 3 , . . ., ZqZ-j and expand any of the bitwistors in 
the numerator as 

(Yj) IJ = a x Z[Zi + a 2 Z{Z( + ... + a^Z] . (2.16) 

The coefficients can be found by contracting with enough bitwistors two get six 
independent equations. More explicitly, one can consider equations of the form 

(Z 2 Z 3 Yj} = a 4 (2345) + a 5 (2356) + a 6 (2367). 



- 14 - 



and solve for the a's. Once this is done, the factor (ABYj) becomes a linear- 
combination of factors in the denominator, thus reducing the degree of the denomi- 
nator and numerator by one. 

The integral in (2.15) is for a general quantum field theory with a planar sector. 
One can continue with the integral procedure in this case but it will take us too far 
away from the main line of the paper. Therefore we concentrate directly on Af = 4 
SYM. In Af = 4 SYM it has been known since the 1990's [26] that all integrals satisfy 
n — m = 4. In modern language, this means that the integrals are dual conformally- 
invariant as discussed in the simple example of the all massless box integral (2.12). 

Iterating the reduction procedure, we can write the any amplitude as a sum over 
pentagons and boxes. But as far as we have seen, the reduction procedure we have 
described so far does not reduce the pentagons any further. Notice that the pentagons 
we have described here are not scalar pentagons, but tensor pentagons — and they are 
manifestly DCI. However, one is always free to choose a basis of bitwistors including 
Y = loo to obtain scalar pentagons, but only at the cost of manifest dual conformal 
invariance. 

But doesn't the reduction procedure of van-Neerven and Vermaseren, when ap- 
plied to Af = 4 SYM, allow for a reduction all the way down to only scalar boxes? 
One might wonder why our analysis so far does not generate this familiar 'box- 
expansion'. The answer is that the reduction to box-integrals is not valid at the level 
of the integrand — only the reduction to boxes and pentagons (scalar or otherwise) 
is valid at the level of the integrand. In order to obtain the ail-too familiar box- 
expansion, it is necessary to parity-symmetrize the integrand — a step that is only 
justified when integrated on a parity- invariant contour, and one which does violence 
to the highly chiral loop-integrands of a quantum field theory such as Af = 4 SYM. 

Here, we should briefly clarify a point which has been unnecessarily confused 
in the literature on Af = 4. Because integrand-level reduction must terminate with 
boxes and pentagons, and box-integrals are both manifestly parity-even and DCI 
while scalar pentagons — which have a factor of (AB 1^) in the numerator — are not 
DCI, the corrections to the box-expansion needed to match the full integrand of 
Af = 4 were first expressed in terms of parity-odd combinations of scalar pentagons. 
This led some researchers to suppose that there was some connection between DCI 
and parity. There is of course no such connection: as evidenced by the extension of 
BCFW to all-loop orders, the full Af = 4 loop-integrand is DCI. 

Especially for theories such as Af = 4 which are DCI, one should strictly avoid 
parity-symmetrization at one-loop or higher. Although scalar pentagon integrals 
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are quite familiar, chiral pentagons are slightly novel — although they have already 
played an important role in the literature (see e.g. [1,27]). The first appearance of 
pentagon integrals occurs for five particles, and there are essentially two possibilities 
that arise: 

2 



(ASF) x (2345>(4512> 



AB 



(AB 12) (AB 23) (AB 34) (AB 45) (AB 51) ' 



where (2 3 4 5) (4 5 12) in the numerator is for normalization 2 and the bitwistor Y is 
simply ZiZ 3 (this is indicated by the dashed-line in the associated figure); and, 



(ABY) xQHSl) 




AB 
5 4 



(AB 12) (AB 23) (AB 34) (AB 45) (AB 51) ' 



where the factor (3 4 51) in the numerator is for normalization, and the bitwistor 
Y = '(512) fl (234)' is the line in twistor-space which lies along the intersection of 
the planes spanned by twistors (Z^, Zi, Z 2 ) and (Z 2 , Z%, Z4) — which is indicated in 
the figure by the 'wavy-line'. As the first of many such examples, it is useful to 
write-out Y explicitly: 

Y = (512) fl (234) = Z 5 Z 1 (223A) + Z 1 Z 2 (5234) + Z 2 Z 5 (12 3 4), 

(2.19 J 

+ Z 1 Z 2 (5 2 34) + Z 2 Z 5 (12 3 4), 

where we have used the fact that (2 2 34) = 0. (The translation between statements 
such as 'the line along the intersection of two planes' and explicit representative 
formulae such as the above will be explained in more detail below; here, we merely 
quote the result in a way from which we hope it will easy to guess the general case.) 

These two integrals are examples of a very important class of integrals that we 
call chiral integrals with unit leading singularities, or pure integrals. In each case, the 
bitwistor appearing in the numerator (together with the integrand's normalization) is 



2 We will see that this normalization follows from the requirement that the integral have unit 
leading-singularities, and its sign is fixed by parity relative to the 'wavy-line' pentagon drawn below 
it. In fact, as we will describe in section 3, the dashed-line in the figure dictates both the bitwistor 
Y = Z1Z3 and the normalization of the integral. 
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completely specified by the dashed- or wavy-line in the corresponding figure. We will 
explain many of the important features of these integrals together with the way their 
graphical representations in more detail in section 3. It is worth noting in passing, 
however, that the two integrals are parity conjugates of one another, and special 
bitwistors Y and Y represent the two lines in twistor-space which simultaneously 
intersect the four lines (51), (12), (23), and (34); this means that (Y 51) = (Y 12) = 
(Y 23) = (Y 34) = 0, and similarly for Y. Because of this, they represent the two 
isolated points in AB-spa.ce for which these four propagators go on-shell. 

Before moving-on to discuss loop integrands, we should emphasize that because 
the primary focus of this paper is the loop integrand — the sum of all the Feynman 
diagrams, as a rational function — there is nothing to say about the regulation of IR- 
divergent integrals such as the zero-mass box integral and the pentagons integrals 
given above. The only integrals we will evaluate explicitly are all manifestly finite 
(in a precise sense which will be described in section 4), and hence are well-defined 
without any regulator. However, it is important to mention that IR divergent inte- 
grals can also easily be regulated and evaluated. In fact, the most natural way to 
add a regulator is also a very physical one, given by moving out on the Coulomb 
branch [28] of the theory. 

2.2 The Loop Integrand 

A simple but far-reaching consequence of writing each Feynman integral in a loop 
amplitude using the dual variables is that one can meaningfully combine all integrals 
appearing in a particular amplitude under the same integral sign. This leads to the 
concept of the loop integrand [1]. We stress again that planarity and the use of dual 
variables plays a crucial role in making this possible-for a general theory, there is no 
natural origin of loop momentum space and therefore no canonical way of combining 
all Feynman diagrams under a common loop integral. 

It is easy to characterize the structure of the n particle 1-loop integrand for Kf = 4 
SYM using momentum-twistor space integrals. All the terms in the integrand can be 
combined defining a universal denominator containing all n physical propagators of 
the form (AB a a+1). If a particular Feynman diagram has fewer propagators, then 
the numerator is chosen so as to cancel the extra propagators. The loop amplitude 
is given as an an integral over a single rational function, 




(2.20) 



AB 
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where A n is the full 1-loop partial amplitude. This formula is already written using 
the simplifications that arise in M = 4 SYM, in other words, it is manifestly DCI. 
However, the integrand exists in any planar theory: for a theory which is not DCI, 
(2.20) would necessarily contain also terms with powers of (AB loo). 

At higher loops, say L loops, scattering amplitudes are given as linear combina- 
tion of integrals of the form 



nti pvt. 



x 



(2.21) 



where N, P, and R are products of Lorentz invariants constructed out of Feynman 
propagators and which depend on the variables shown and on the external momenta. 
Written in this form, there is clearly a large amount of redundancy in the definitions 
of the internal loop momenta. 

Since we are dealing with only planar integrals, for each Feynman diagram there 
exists a dual diagram (the standard dual graph of a planar graph). Consider for 
example the following four-point two-loop integral: 
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(2.22) 



Using Xi to denote the dual coordinates of the external momenta and yi to denote 
the internal points, one can write any planar L-loop integral in dual coordinates. 
There is, however, one slight subtlety in using such a prescription to uniquely define 
'the' integrand: while the definition of the external points is unique, the labeling 
of the internal points is not (when L > 1). But the solution to this problem is 
very simple: we are always free to completely symmetrize the integrand with respect 
to all L\ permutations of the internal loop- variable labels. Although we will often 
write multi-loop integrands in some particular representative choice of the labels for 
internal propagators, complete-symmetrization over all permutations of indices is 
always implied (including a factor of 1/L! from this symmetrization). 

Consider for example the simplest two-loop integral, given above in (2.22). Writ- 
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ten in dual-coordinates, the integral would be given by 

d 4 yid%\ ((xi~ x 3 ) 2 ) 2 (x 2 - x 4 ) 2 



+ (yi 



2 J (Vi- ^3) 2 (yi - x 4 ) 2 (yi - x 1 ) 2 (y 2 - x 1 ) 2 (y 2 - x 2 ) 2 (y 2 - x 3 ) 2 (y 1 - y 2 ) 2 

—where the numerator was chosen in order to make the integral dual-conformally 
invariant, and the factor of 1/2 in the measure reflects the complete-symmetrization. 

Of course, as we will see repeatedly throughout this paper, (multi-)loop inte- 
grands are much more naturally expressed in terms of momentum-twistor variables. 
To translate the integral (2.22) in momentum-twistor variables, we need to associate 
a pair of twistors to each of the two loop variables. This we can do by making the 
association 

yi ^(Z A ,Z B ) and y 2 ^(Z c ,Z D ). (2.23) 

Using this notation and the translation of propagators in terms of momentum 
twistors given in (2.5) one finds 

(1234) 2 (2341) 



(AB 41) (AB 12) (AB 23) (CD 23) (CD 34) (CD 41) (AB CD) ' 

(AB,CD) 

where '(AB, CD)' implies that the integration measure carries with it a factor of 
1/2 from the symmetrization of (AB) -H- (CD). We should mention here that for 
3-loops, we will use (Ze, Zp) to denote the line corresponding to y 3 — but of course, 
a convention such as that of associating (Z \ m ,Z Bm ) with y m would be increasingly 
preferable at high-loop order. 

Before we leave the topic of the loop- integrand in general, we should mention that 
the form of the integrand obtained via BCFW as described in [1] makes it completely 
manifest that the loop-integrands in M = 4 enjoy the full Yangian symmetry of 
the theory. (Of course, the choice of an integration contour which introduces IR 
divergences, such as the physical contour, breaks this symmetry.) 

However, just as with the BCFW recursion relations at tree level, the formulae 
obtained from the recursion do not enjoy manifest locality or manifest cylcic invari- 
ance. The restriction that we impose throughout this work, however, is that loop- 
integrand be expanded in a way which makes use of only planar, local propagators. 
As we have stressed a number of times, we will find amazingly simple, manifestly 
cyclically symmetric and local expressions for multi-loop amplitudes, that are signif- 
icantly simpler and more beautiful than their BCFW counterparts! Taken together 
with the parallel results presented in [9], this strongly suggests the existence of a 
formulation for scattering amplitudes directly yielding these remarkable local forms. 
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The local formulae presented in this paper are very closely related to and influ- 
enced by the concept of the leading singularities of scattering amplitudes, which we 
proceed to presently describe. 



2.3 Leading Singularities 
Definition 

The concept of leading singularities was introduced in the 1960's in the context of 
massive scalar theories [29]. More recently, in 2004, the same concept was modified 
to accommodate massless particles and this was exploited for Yang-Mills in [30]. 
The original definition of 'leading-singularity' refers to a discontinuity of a scattering 
amplitude across a singularity of the highest possible co-dimension. At 1-loop, for 
example, leading singularity discontinuities are computed using a generalization of a 
unitarity cut, but where four propagators are cut instead of two. Using A% for i = 
1, . . . , 4 to denote the four partial amplitudes, each with their associated momentum- 
conserving 5-function, one has what can be called leading-singularity discontinuity, 




(2.24) 



Y[d%d%5{£ 2 r ) A 1 {{£i,m}A-t2^>---)xM{t2MA-tsM, 



r=l 



x A 3 ({£ 3 ,V3},{-^v4,---) x 774}, {-^i,77i}, •••)■ 



Here, the integrations over the internal loop momenta are there only to remind us 
that we are to sum-over all solutions to the conditions imposed by the (^-functions, 
and the integral over the Grassmann coordinate rji of each internal particle £j is 
there to remind us that we are to sum-over the exchange of all possible internal 
particles — which in the case of M = 4 means the full super-multiplet. 3 



3 Here, we are using an on-shell superspace formalism which allows us to talk about all particles 
in the same super-multiplet as a single 1-particle state. We assume familiarity with this concept, 
but for careful definitions, more references and applications see [31]. 
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This point of view of leading-singularities has been very useful and allows a 
complete determination of 1-loop amplitudes in M = 4 and in M = 8 supergrav- 
ity amplitudes when thought of as linear combinations of scalar box integrals with 
rational coefficients. The rational coefficients can be computed using the notion of 
generalized unitarity. Clearly, the notion of discontinuities is not related to the exis- 
tence of an integrand and this is the reason it works in Af = 8, supergravity where an 
analog of 'the integrand' — which requires a way to combine integrals with different 
cyclic orderings — has not yet been found. 

As mentioned in our discussion of reduction procedures in M = 4 SYM, the 
expansion in terms of boxes cannot give the physical integrand. The physical inte- 
grand is defined as that which coincides with the one from Feynman diagrams, prior 
any to reduction techniques, as rational functions — and, as we will see, the Feynman 
diagrams of M = 4 in a given i?-charge sector are chiral. 

Once we think about the integrand as being the object we are after, we can 
try to model it by using some appropriate basis of functions, dictated by a general 
reduction procedure. Clearly, the set of all DCI tensor pentagons and boxes should 
be enough. Nevertheless, we will find that such a basis would still possess many of 
the unattractive features of the box-expansion, and so we will introduce much more 
refined choice in section 5. 

The importance of dealing with a specific rational function is that we can inte- 
grate it on any choice of contour we'd like — not just the real-contour which defines 
the Feynman integral. This allows us to define a more refined notion of a leading- 
singularity — the previous notion, motivated by generalized unitarity, is much coarser 
version of the one we will use now. In [32], this more refined notion was introduced, 
and it was used to match the full M = 4 integrand for several 1-loop and 2-loop 
examples. However, in [32] the deep reason for why the idea was working, i.e., the 
existence of the integrand, was not appreciated. 

Whether written in ordinary momentum space, using dual-coordinates, or using 
momentum-twistors, loop integrals can be thought of as complex contour integrals 
on C 4 with the choice of contour corresponding to R 4 — the real-slice. However, this 
choice of contour is known to break many of the symmetries of the theory, and is 
littered with IR-divergences, etc. that can be the source of confusion. From various 
viewpoints, the most natural contours would instead be those which compute the 
residues of the integrand. These are always finite, are often vanishing, and make 
manifest the full Yangian symmetry of the theory. We refer the reader to [33] for 
a mathematical definition of residues in several complex variables; here we hope 



- 21 - 



the reader will find the definitions a natural generalization of the one-dimensional 
residues with which everyone is familiar. 

Let us present the definition using x variables first. Consider a contour of inte- 
grations with the topology of a T 4 = (S* 1 ) 4 . In order to compute a particular residue 
one has to choose four propagators (x — x ai ) 2 , with i = 1, ... ,4 and integrate over 
the T 4 , defined by \{x — x ai )\ = where are small positive real numbers near 
one of the solutions. The circles, S 1 are parametrized by the phases and are given a 
particular orientation. 

The definition of a multidimensional residue is very natural if one defines vari- 
ables Ui = (x — x ai ) 2 - Performing the change of variables the integral becomes 



where now the contour becomes small circles around Ui = 0. J is the Jacobian of the 
change of variables. The residue is then the Jacobian times the rest of the integrand 
evaluated at Uj = 0. The Jacobian 



is clearly antisymmetric in the order of the columns. Different orderings can differ by 
a sign and this is related to the orientation of the contour. These signs are important 
when discussing the generalization of residue theorems to the multidimensional case, 
which will play an important role momentarily. 

From now on we call each individual residue a leading-singularity. As before, 
these are given by the product of four on-shell tree amplitudes as shown in Figure 2. 




x — x {The rest of the integrand} 



(2.25) 




(2.26) 




Figure 2: A 'quad- cut' one-loop leading-singularity viewed as a T 4 contour-integral which 
'encircles' the point in C 4 where four-propagators are made on-shell. 



The reason for the appearance of the tree amplitudes is that the residue of the poles 
is computed where the four propagators vanish and therefore internal particles can 
be taken on-shell. 

Leading singularities at higher loop-level can also be defined as residues of a 
complex, multidimensional integral over C 4L where L is the loop order. This means 
that in order to define a residue one has to define a T 4L torus as a contour of integra- 
tion. Naively, residues can only be defined for integrals with at least 4L propagators. 
However, noticing that propagators are quadratic in the loop-momentum, one can 
define composite leading singularities which involve less than 4L propagators as done 
in [4,32,34], using the self-intersection of curves defined by the on-shell condition to 
define isolated points in C 4 about which the T 4i contour should 'encircle.' 

We will not discuss composite leading singularities in detail here simply be- 
cause we will present evidence that when a special set of integrals, we call chiral 
integrals with unit leading- singularities, are used, matching non-composite leading- 
singularities appears to suffice to fix the entire amplitude. Moreover, we will see that 
only a very small subset of non-composite leading-singularities need to be considered 
to accomplish this. 

Chirality of Leading Singularities 

It turns out that for nonsingular external momenta, there are exactly two solutions 
to the equations (x — x a .) 2 = 0, with % — 1, . . . , 4, and therefore two residues of each 
choice of four propagators. (This has a beautiful geometric interpretation in momen- 
tum twistors as we will see shortly.) This means that for an n-particle amplitude, 

there are 2^ (non-composite) one-loop leading-singularities. 

Consider any box integral, say, an integral with two massless legs and two mas- 
sive, known as the 'two-mass-easy' integral: 
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(2.27) 
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where N is just some normalization that need not concern us presently. The equa- 
tions 

(x — Xi_i) 2 — (x — Xi) 2 — (x — Xj_i) 2 — (x — Xj) 2 = 



have two solutions, and therefore a residue can be computed for each such point 
separately. We'll soon see that these two solutions are easily found and differentiated 
when written with momentum-twistor variables; but for now, let us suppose the two 
solutions have been found, and denote the corresponding contours Tf and T 2 4 . 

A very important tool that will make an appearance many times is multidimen- 
sional analogue of Cauchy's theorem, called the Global Residue Theorem (GRT). 
The GRT states that — given a suitable condition at infinity — the sum over all the 
residues of a given rational function vanishes (see chapter 6 of [33]). This means, in 
the present case, that 

res T 4(X 2me ) + res T 4(X 2mc ) = (2.28) 

Moreover, we can choose the normalization N is such that, say res T 4(X 2mc ) = 1. Such 
a choice is possible for all box integrals, following from the simple fact that all box- 
integrals — having only four propagators — must have residues which are proportional 
equal and opposite. We refer to this fact by saying that scalar box integrals are not 
chiral. The use of the word chiral is justified by the fact that the locations of the 
leading singularities, as points in C 4 , are mapped into each other by parity — which 
is just complex conjugation. And so the corresponding contours are mapped into 
each other up to orientation by parity. If use (T 4 )* to denote the parity conjugate 
contour of T 4 , then res,- T 4), = — res T 4 and the GRT implies that 

res T 4(X 2mc ) = res (T 4)*(X 2mc ). (2.29) 

Let us now consider the leading-singularities of the one-loop integrands of M = 4 
Yang-Mills. We'll see that, as scattering amplitudes of Af = 4 in a given i?-charge 
sector are chiral, so are the one-loop leading-singularities of field theory! In other 
words, the two residues associated with the two solutions of cutting four-propagators 
are not the same. Let us see this in an example. The simplest possible example is 
the five-particle MHV amplitude 4 . Let us consider taking the leading singularities 
of the field-theory integrand which encircles the point in C 4 where the following four 



4 The only DCI object for four-particles is the zero-mass box integral. This is why both leading 
singularities are equal to the tree amplitude. 
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propagators go on-shell: 




(x — xi) 2 — (x — X2) 2 — (x — X3) 2 — (x — X4) 2 = 0. (2.30) 



It was noticed already in [30] that on one solution M = 4 SYM gives the tree 
amplitude, A^ ee , while it vanishes on the second. 

The vanishing of leading singularities can be understood from pure supersym- 
metry. Consider an amplitude in the i?-charge sector m. Recall the N m ~ 2 MHV 
classification of amplitudes in J\f = 4: under a rescaling of all rj a variables by trf a , 
an N m_2 MHV amplitude picks up a factor of t 4m . From the definition of leading 
singularities as the product of tree amplitudes connected by internal on-shell states 
we see that every internal line contributes (—1) to the .R-charge counting coming 
from the integration over rj variables. At 1-loop, we have four tree-amplitudes and 
four propagators. If the -R-charge of each tree-amplitude is rrii (see Figure 2), then 
the -R-charge of the leading singularity is m\ + + + — 4. 

Returning to the five-particle example, because we are interested in a one-loop 
MHV amplitude, all its leading-singularities must have m = 2. The four-particle 
vertex (in the upper-left of the figure above) can only have mi = 2 and therefore the 
three-particle vertices have to satisfy mi + 1713 + = 4. Since the possible values 
for m for a three-particle amplitude are 1 and 2, two vertices must have m = 1 and 
one must have m = 2. This leaves only the possibilities shown below: 



2 2 : 



Of these three possible leading-singularities of field theory, it turns out that the first 
one is equal to the five-point MHV tree-amplitude, and the latter two vanish for 
generic external momenta. In fact, whenever one is considering a leading singularity 
which involves 3-particle vertices, some very simple and powerful rules prove very 
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useful: 1. any leading singularity involving adjacent three-particle vertices with the 
same i?-charge will vanish for generic external momenta (momentum conservation 
in this case, requires that the external particles attached to these vertices must be 
collinear); and 2. leading singularities involving three-particle amplitudes are almost 
always chiral — the only exception being the four-particle amplitude. 

In the case of the five particle example under consideration, we see that the 
residue from the contour encircling one of the two solutions to the quad-cut equations 
in (2.30) is equal to *4.5mhv> wn ^ e the conjugate contour integral vanishes. We will 
explore this in more detail once we introduce the geometric point of view. 

Dual Formulation of Leading Singularities 

In the rest of the paper we will make much use of the fact that leading-singularities 
satisfy many relations. These relations can be seen as resulting from residue theorems 
of the integrals which compute them. As a final comment before exploring the 
connection between leading singularities and the classic enumerative problems in 
the projective geometry of momentum twistor space let us briefly introduce the 
Grassmannian formulation. 

In [4], leading singularities were proposed as completely IR-finite quantities that 
were likely to contain all the information needed to compute the S-Matrix of M = 4 
SYM. Moreover, it was conjectured that all leading singularities of the theory, which 
can be obtained to arbitrarily higher loop order, are computed by a contour integral 
over a Grassmannian manifold 5 G(m, n) called C m , n - Here m determines the R- 
charge sector of the theory under consideration. 

The integral was first presented in twistor space 



In this presentation, residues of this integral are manifestly superconformal invariant 
(that is, superconformally-invariant in ordinary spacetime). Here we have introduced 
the concept of dual super twistor space W = (A, /i, rj). This particular space will not 
play a significant role in this work, so we refer the interested reader to [4,35] for more 
details. 

This formula can be transformed to momentum-space and then to momentum- 
twistor space. Very remarkably, the formula in momentum-twist or space also turns 

5 The Grassmannian G(m, n), a natural generalization of ordinary projective space, is the space 
of m-dimensional planes in n-dimensions. Each point in G(m,n) can be represented by the m 
n- vectors which span the plane, modulo a GL(m) redundancy. 





(2.31) 



vol(GL(m)) (12 ••• m)(2 3 ••• m+l)---(nl ••• m-1) 
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out to be an integral over a Grassmannian, with the MHV-tree-amplitude arising as 
the Jacobian from the change of variables. Specifically, 

£m,n Imomentum— space (^j A ( Tj) C-2,n X T^k,n: (2.32) 

where k = m — 2 and 

mW ~ 7 vol(GL(fc)) (1 2 • • • k){2 3 • • • fc+1) • • • (n 1 • • • fc-1) " 1 J 

This representation in momentum twistor space makes dual superconformal invari- 
ance manifest [36,37]. With some more effort one can prove that residues of this 
formula are also invariant under level one generators of the Yangian of the dual su- 
perconformal algebra and hence invariant under the whole Yangian [38]. The level 
one generators are nothing but the superconformal generators when passed through 

It has now been proven that all leading singularities are Yangian invariant and 
that all Yangian invariants are residues of the integral (2.33). From the physical 
point of view the problem has been solved. It might also be interesting to go further 
and prove that all residues of (2.33) correspond to some leading singularity but we 
will not discuss this issue any further. 

Momentum Twistors and Schubert Problems 

Statements like the number of solutions to setting four propagators to zero is two are 
non-obvious from the dual space x point of view. In terms of momentum twistors, 
this statement turns out be a simple, classic problem of the enumerative geometry 
of CP 3 , solved by Schubert in the 1870's [39,40]. 

Recall that an n-particle 1-loop amplitude can be written as 

A f Y, i c i (ABY;)(ABY})---(ABYi_ A ) 
n J (AB 12) (AB 23) ■■■(ABn-ln)(ABnl)' 1 j 

AB 

Each one-loop leading-singularity is associated with a point in the space of loop- 
momenta for which some choice of four propagators simultaneously become on-shell, 

i i+l 

cL cb (ABi i+l) = (ABj j+1) = (AB k = (ABl Z+l) = 0; 

fc 
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Because the loop momentum is represented in momentum-twistors as the line (AB), 
the solution to these four equations should correspond to a particular configuration 
for the line (AB). We will see that for all leading-singularities which involve a three- 
particle vertex (a 'massless leg'), the two solutions to four equations above are cleanly 
distinguished geometrically, allowing for a richly-chiral description of the integrand. 

Before describing the full problem of putting four propagators on-shell, let us 
briefly consider the geometric significance of having a single factor, say (ABii+1), 
vanish. Recall that the four-bracket (••••) is nothing but the determinant of the 4x4 
matrix of components of its four momentum-twistor arguments (viewed as elements 
of C 4 ). As such, (ABii+1) = if and only if the vectors Za, Zb, Z{, Zj + i are not 
linearly independent, implying the existence of some linear relation among the four 
twistors of the form (XaZa + olbZb + CKjZj + a i+ iZ i+ i = 0. Trivially rearranging we 
see that 

c*aZa + a B Z B = -(a.iZi + a i+ iZ i+1 ), (2.35) 

which we may read as saying there is a point on the line spanned by Za, Zb — namely 
(c^aZa + ol b Zb) — which lies along the line spanned by Zj, Z i+1 . Which is to say, the 
lines (AB) and (ZjZ i+1 ) intersect; and because two intersecting lines describe a plane, 
we say that the four points Za, Zb, Z^ Z i+ \ are coplanar. 

Therefore, the problem of finding the particular lines (AB) for which four prop- 
agators simultaneously vanish is equivalent to finding the set of lines in CP 3 which 
simultaneously intersect four given lines (which are presumed fixed by the external 
data). The number of solutions to this problem is one of the classic examples of the 
enumerative geometry developed by Schubert in the 1870's. For this reason we call 
these problems Schubert problems. 

The answer to the number of lines which intersect a given four turns out to be 
remarkably robust: provided the four lines are sufficiently generic, there are always 
2 solutions, and an infinite number otherwise. 6 (An example of a non-generic con- 
figuration would be one for which three or more of the lines were coplanar; these are 
never found for generic external momenta.) 

Schubert derived the number of such solutions with an argument that is decep- 
tively simple. The idea is to consider a particular configuration where it is easy to 
count the number of solutions. Schubert intuited that the answers to such enumer- 
ative questions should be topological in nature, and therefore should not depend 
on the particular configuration in question. Therefore, one can analyze the most 

6 To be precise, we must count solutions with multiplicity; however, for a generic set of lines in 
the problem, the 2 solutions will always be distinct. 
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convenient possible configuration (for which the number of solutions is not infinite) 
and the answer found for that case, should be the answer in general. Said another 
way, the number of solutions to a given Schubert problem should not change when 
a particular special configuration is smoothly moved into a more general position. 

Perhaps the easiest configuration for which we can count the number of solutions 
to the Schubert problem of finding the lines (AB) that intersect four given lines in 
CP 3 is the zero-mass configuration; it is so-called because it is the configuration 
which corresponds to the box integral with zero of its four corners massive, 



(1234) (2341) 
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1 4 

which is an integral we have seen before. Explicitly, we would like to find all the 
lines (AB) which intersect all the four lines (12), (23), (34), and (41). This problem 
is indeed easy to solve, and the two solutions are drawn below. 





2 

(AB) = (13) 

Clearly, because (12) fl (23) D Z 2 and (34) fl (41) D Z 4 , the line (AB) = (24) inter- 
sects all four lines, as desired; this is drawn in red above. The same argument also 
applies to the second solution, the line (AB) = (13), drawn in blue above. Also in 
this figure, we have indicated which leading-singularities have non-vanishing support 
on the corresponding (complex) point in the space of loop-momenta which corre- 
sponds to the particular line (AB). As explained above, each three-particle MHV 
(m = 2) — colored blue in the figure above — or MHV (m = 1) — colored white — vertex 
of a leading singularity vanishes for every leading-singularity, and so which of the 
2 three-particle amplitudes is non-vanishing for this value of the loop-momentum 
determines the chirality of the contour. 

As a convenient way to gain some intuition about momentum-twistor geometry 
that will prove useful in the rest of this paper and to establish some of the notation 
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that will be ubiquitous throughout, we will study each of the 1-loop Schubert prob- 
lems in turn. 



One-Mass Schubert Problem: 

A 'one-mass' 1-loop leading singularity is one for which three of the four legs are 
massless, and is associated with the following archetypical box-integral: 
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W 
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AB 



(AB 12) (AB 23) (AB 34) (AB 45) 



(2.36) 



In momentum-twist or space, the leading-singularities of this integral are associated 
with the lines (AB) which intersect the four lines (12), (23), (34), and (45). Consid- 
ering the configuration of lines, it is not hard to find the two configurations which 
solve this Schubert problem: 



2 3 
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(AB) = (24) 




(AB) = (123) H (345) 



As before, because (12) fl (23) D Z 2 and (34)0(45) D Z 4 , the line (AB) = (24) 
intersects all four lines. The second solution, however, is new. This solution is drawn 
in blue in the figure above, and represents the line of the intersection of the planes 
spanned by (Z%, Z 2 , Z 3 ) = (123) and (Z 3 , Z 4 , Z 5 ) = (345). Although geometrically 
clear, it is worthwhile to recall that any generic line in the plane (123) will intersect 
the lines (12), (23), and (31), and any generic line in the plane (345) will intersect 
the lines (34), (45), and (53). Therefore, the line (AB) = (123) fl (345) will intersect 
all four lines, as required. 

Similar to the case discussed in the context of the pentagon with a 'wavy-line' 
numerator (2.18), the line (123) fl (345) can easily be expanded in terms of ordinary 
bitwistors as: (23) (1 345) + (31) (2 345). This follows from a more general rule which 
review presently. 
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On the Intersection of Planes in Twistor- Space 

In general, the intersection of the planes (abc) f] (def) is can be canonically 
expanded in either of the following ways: 

(abc) (1 (def) =Z a Z b (cdef) + Z b Z c (adef) + Z c Z a (b def); 

= {abcd)Z e Zf + (abc f)ZdZ e + (abce)ZfZ ( i- 

Alternatively, when expanding a four-bracket of the form (xy (abc) f] (def)), the man- 
ifest dependence on the two planes can be preserved at the cost of breaking the 
manifest dependence on the line (xy), as follows: 

(xy (abc)f](def)) = (x abc) (y def) — (y abc) (x def). (2.38) 



Two-Mass-Easy Schubert Problem 

The two-mass-easy Schubert problem is associated with the following one-loop 
archetypical box-integral, 
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(2.39) 



which has leading singularities supported on the configuration (AB) which intersect 
all four of the lines (12), (23), (45), and (56). The two solutions are essentially the 
same as for the one-mass Schubert problem, and are illustrated in the Figure below: 
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(AB) = (25) 



(AB) — (123) H (456) 



Once again, there is a very easy solution, in this case the line (AB) = (25) which 
obviously intersects the four lines. And using the same reasoning as int the one-mass 
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Schubert problem, it is easy to see that the second solution is simply the intersection 
of the planes (123)0(456). 



Two-Mass-Hard Schubert Problem 

The two-mass-hard Schubert problem differs from the two-mass easy problem in 
that the two massless corners are adjacent — making the Schubert problem slightly 
less 'easy' (which at least partially justifies the name). It is associated with the 
following archetypical one-loop integral, 
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(2.40) 



and has leading singularities supported where the line (AB) intersects the four lines 
(12), (23), (34), and (56). The two solutions are shown in the Figure below: 




Let us briefly discuss the first of the two solutions. Here, the line (AB) = (123) f) (356) 
intersects the lines (23), (34) trivially because Z 3 C (123)0(356), and it intersects 
the lines (12) and (56) because any generic line in the plane (123) intersects (12), 
and any generic line in the plane (356) intersects (56). 



Three-Mass Schubert Problem 

The last Schubert problem that involves a massless corner is known as the 'three- 
mass' problem, and is associated with the following archetypical one-loop integral: 
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This integral is the most general one which involves a massless corner, and supports 
leadings singularities where the line (AB) intersects the four lines (12), (23), (45), 
and (67). The two solutions are indicated in the Figure below. 
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(AB)= ("(123)0 (45), (67)0(123) 
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(AS) = (245) (672) 



Here, the notation '(ab) f] (cde)' has been used to indicate the point in twistor-space 
where the line (ab) intersects the plane (cde). We will discuss the expansion of such 
geometrically-defined objects more generally at the end of this subsection; for now, 
let us merely quote the result: 

(ab) P| (cde) = Z a (bcde) + Zi(cdea) = —(^Z c (deab) + Z c i(eabc) + Z e (abcd)j] 

and similarly, 

(cde) Pi (ab) = Z c (dea b) + Zd(e a be) + Z e (abcd) = — (^Z a (bcde) + Z^cdea)^ ; 



so that (ab) f] (cde) = —(cde) f] (ab). 



On Schouten-Identities and Projective Geometry 

Perhaps the single most useful identity for momentum-twistor geometry is known 
as 'the five-term identity:' any arbitrary set of five twistors {Z a , Zf,, Z c , Z^, Z e } will 
satisfy the following identity, 

Z a (bcde) + Z b (cdea) + Z c (deab) + Z d (eabc) + Z e (abcd) = 0. (2.42) 
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This identity merely reflects the general solution to a homogeneous, linear system of 
equations in four-variables, and as such, has analogues in any number of dimensions. 
For example, in two dimensions, we have that for any {A a , Ab, A c } C C 2 , there is an 
identity 



where we have naturally extended the definition of '(■•)' to be the determinant of 
the components of the corresponding two-vectors. This two-dimensional identity 
represents the general solution to a homogeneous, linear system of equations in 2 
unknowns, and by contracting it with a fourth two-vector A^, we obtain the familiar 
'Schouten identity:' 



This familiar identity of course has an analogue descending from equation (2.42). 
By contracting equation (2.42) with any arbitrary plane (/ g h), we find the following 
5-term identity which we will therefore call 'a Schouten identity:' 

(/ g ha)(bcde) + (f g hb)(cdea) + (f g hc)(deab) + (f g hd)(eabc) + {f g he) (a bed) = 0. 

In addition to being quite useful for simplifying formulae, equation (2.42) can 
be trivially re-arranged to yield the solutions to some of the most often-encountered 
problems in momentum-twistor geometry: 

1. the expansion of any arbitrary twistor Z a into a basis composed of any four 
linearly-independent twistors {Z b , Z c , Z d , Z e }: 



^ a {bc) + Afc(ca) + \ c (ab) = 0, 



(2.43) 



(da){bc) + (db)(ca) + (dc){ab) = 0. 



(2.44) 




Z b (cdea) + Z c (deab) + Z d (eabc) + Z e (abcd) 



2 



the point along the line (ab) which intersects the plane (cde): 





3 



the point on the plane (abc) which intersects the line (de): 



(abc) fl (de) = Z a (b c d e) + Z b (c d e a) + Z c {d e ab) 




and so-on. 
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Matching All Leading Singularities 

We close this introductory section to momentum twistor integrals and leading singu- 
larities with a physical point. We have seen that the leading singularities of M = 4 
SYM are chiral while those of scalar boxes are non-chiral. This means that if we 
want to construct the integrand of the theory it is impossible to do it using scalar 
boxes. Momentum twistors already give the solution to this problem. Since leading 
singularities are Yangian invariant and in particular dual conformal invariant (DCI), 
one should use the reduction procedure to go down to tensor pentagons and boxes 
and not any further. Even going down to scalar pentagons would be doing something 
brutal to the manifestly DCI structure of the amplitudes. 

In the rest of the paper we will find that by using a special class of integrals 
known as chiral unit leading singularity integrals, the full integrand of scattering 
amplitudes can be reproduced yielding to stunningly simple forms. 
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3. Chiral Integrals with Unit Leading Singularities 

Given the success of the recently introduced recursion relations for the construction 
of the integrand to all orders in perturbation theory [1], it is clear that the physical 
integrand is the important object to obtain. 

In the previous section we showed that the usual constructions of, say, one-loop 
amplitudes in M = 4 SYM as a linear combination of scalar boxes cannot possibly 
be the physical integrand. Of course, the answer obtained from scalar boxes gives 
the same integrals as the one originally defined from Feynman diagrams. However, 
as we will see, insisting in obtaining the physical integral leads to stunningly simple 
formulas for one and higher loop amplitudes. These new formulas are possible thanks 
to the use of a new suit of integrals with very special properties. These are chiral 
integrals with unit leading singularities. 

3.1 Integrals with Unit Leading Singularities, or Pure Integrals 

Let us start by given a definition of integrals with unit leading singularities. As we 
will see, it is appropriate to call these pure integrals. 

Consider a particular DCI L-loop integral and compute all possible residues. 
If all non-vanishing residues are the same up to a sign then the integral can be 
normalized so that all residues are ±1 or 0. When this is done, the integral is said 
to have unit leading singularities or to be a pure integral. 

We already encountered examples of pure integrals in the previous section. The 
zero mass box (2.12), the general scalar box (2.13) (properly normalized), and the 
pentagon integrals in (2.17) and (2.18). 

Using the global residue theorem, we proved in section 2 that boxes are pure 
integrals. However, it is not obvious that the pentagons in (2.17) and (2.18) satisfy 
the requirement. 

Consider first pentagons of the first class 



where N = (12 45) (23 45). 

In order to see that all non-vanishing leading singularities are equal up to a sign 
let us use a global residue theorem. In section 2 we gave a very imprecise definition 
of the global residue theorem (GRT) which was enough for the purposes of that 
section. Here we have to be more precise. The GRT states that given a choice of a 
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map / : C 4 — > C 4 made from polynomial factors in the denominator, the sum over 
all the residues associated with the zeroes of the map vanishes. 

In the present case, consider the map given by / = / 2 , f'3, f&) where 

h = (AB 12), h = (AB 23), f 3 = (ABM), U = (AB45)(AB51). 

It is easy to see that the map / has four zeroes (see section 2 for more details). The 
GRT assures that the sum over the four residues vanishes. How can we prove that 
residues are equal if the GRT only gives relations among four residues? 

The answer has to do with our choice of numerator. Consider the value of 
(AB 13) on the four zeroes. Each zero is a line which is the solution to some Schubert 
problem 7 . The four solutions are the lines (24), (123) fl (345), (13) and (512) fl (234) 
(see the end of the section or section 2 for the notation). It is a simple exercise to 
show that (AB13) vanishes on the second and third solutions and it is non zero on 
the first and fourth. This means that the GRT implies that two leading singularities 
are equal and opposite in sign. The first is one of the two solutions to (AB 12) = 
(AB 23) = (AB 34) = (AB 45) = while the fourth is one of the two solutions to 
(AB 12) = (AB23) = (AB 34) = (AB 51) = 0. Let us denote these non-vanishing 
residues by r ( i 2 ),(23),(34),(45) and r (12 ) i ( 2 3),(34),(5i) respectively. Therefore the GRT states 
that 

(0 + r(i2),(23),(34),(45)) + (?"(12), (23), (34), (51) + 0) = 

which implies the equality of the residues up a sign. 

The pentagon integral as 10 leading singularities. This means that more work is 
needed to show that it has unit leading singularity. Consider a GRT associated to 
the map 

h = (AB12)(AB51), h = (AB23), f 3 = (AB34), / 4 = (AB45). 

Once again, there are four zeroes of this map. Two of them are shared with the map 
we constructed before, i.e., (24) and (123) H (345). The two new solutions are (35) 
and (234) f] (451). As before, the numerator vanishes on (123) f] (345). Very nicely, 
it also vanishes on (35). We can denote by r(i2),(23),(34),(45) and r (51), ( 2 3), (34), (45) the 
corresponding non-zero residues. Therefore the GRT gives 

(0 + r(i2),(23),(34),(45)) + (j(51),(23),(34),(45) + °) = 

7 A Schubert problem was defined in section 2 as the projective geometry problem of finding 
lines that intersect four given lines which can be in special configurations called one-mass, two- 
mass-easy, two-mass-hard, and three-mass, as well as in generic positions which we call four-mass 
configurations. 
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This means that the GRT sets equal the non vanishing leading singularity in (AB 51) = 
(AB 23) = (AB 34) = (AB 45) = with the ones we found before. 

This procedure can be continued three more times by shifting the labels in the 
map by one. We leave it as an exercise for the reader to verify that in every case, the 
numerator vanishes on one solution implying that the GRT sets all non-zero leading 
singularities to be the same. 

In order to compute the normalization and also to show how the GRT makes 
obvious statements that require computations to be verified, even in this trivial case, 
let us compute explicitly the two residues in the first GRT discussed above. 

Consider the ones in the first step. In other words, let's evaluate the residue on 
the solution (24) to the system (AB 12) = (AB 23) = (ABM) = (AB 45) = 0. The 
residue is given by 

N « (3 2) 

(2451) ((1234) (2345)) 1 ; 

Here the terms in parenthesis are the Jacobian in the residue computation. A ge- 
ometric way to see that the Jacobian has to contain the factors (1234) and (2345) 
is that on the special configurations where either one of them vanishes, the num- 
ber of solutions to the Schubert problem becomes infinite. For example, consider 
the configuration where (1234) = 0. In this case, any line on the plane (123) which 
passes through Z 4 solves the Schubert problem. Using the scaling of each momentum 
twistor, the Jacobian must be what we found. It might be instructive to see the full 
computation of the jacobian using momentum twistors. This is carried out in detail 
in appendix A. 

In order to have a properly normalized integral we require (3.2) to be equal to 
one. This means that A" = (5124) (2345) which is the factor first given in section 2 
in (2.17). 

Consider now the residue coming the second Schubert problem, (AB 12) = 
(AB 23) = (AB 34) = (AB 51) =0. The non-zero residue is associated with the solu- 
tion (512) Pi (234). This is a one-mass Schubert problem and one explicit form of Za 
and Z B was given in section 2. Let us use Z A = Z 2 and Z B = — (1234) Z 5 + (5234) Z\ 
and compute the residue. The Jacobian is the same as before but with labels shifted 
back by one. The residue is then 

N (1234) (2513) (3 3) 

(2345)(5124)((1234)(2513))' { ' 

Using the normalization derived above this quantity equals one as expected. 
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In section 2 we also presented a second pentagon integral which differs from the 
first one only in the choice of numerator. We leave it as an exercise for the reader to 
repeat the analysis done here and show that with the new numerator this is a pure 
integral 8 . Let us rewrite the integral here with the numerator given in geometric 
form 



Now it should be obvious that the comment made in section 2 is true. The special 
numerators are made from lines, (13) and (512) f] (234), which are the two solutions 
to a Schubert problem. 

In section 4 we study a less trivial example; a hexagon integral where the special 
choice of numerator also allows the use of the GRT to show that all non-vanishing 
residues are equal. In the hexagon case, checking the statement that all residues are 
equal algebraically requires many applications of 4-bracket Schouten identities. 

Basic Diagrammatic Notation 

We find it convenient to introduce a diagrammatic representation for numerators. 
Note that with our definition of dual variables p a = x a — x a _i and of momentum 
twistors x a -H- (Z a ,Z a+ i), there is a natural diagrammatic relation between loop 
integrals and momentum twistor configurations. Consider a general one-loop ampli- 
tude as a polygon with n-sides. Attached to each vertex there is some momentum 
p a . In momentum twistor space, we also have an n-sided polygon and attached to 
each vertex there is a momentum twisor Z a . Following the intuitive correspondence 
between the two diagrams we are led to denote denominators (propagators) as lines 
connecting points depending on their geometric configuration. These are denoted by 
solid lines. In order to distinguish numerators, we also introduce dashed and wavy 
lines. 

Dashed lines: Numerators which correspond to factors of the form (ABef), 
where (e/) represents a line in momentum twistor space specified by two momentum 

8 Of course, one could simply translate the whole problem into dual momentum twistor space 
to find exactly the same integral as before. However, it is still an instructive exercise to do it in 
momentum twistor space. 




(3.4) 
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twistors Z e and Zf is represented by a dashed line connecting points e and / as in 
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(AB 13) (1245) (2345) 



: / (Afil2)(Afi23)(Afi34)(Afi45)(Afi51) 
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I^avy /mes: We also allow points to represent dual twistors. In this case the 
second class of numerators constructed as intersection of planes can also be repre- 
sented by a line connecting two points. In order to distinguish this from the previous 
case we use wavy lines. In the example where the numerator corresponds to the line 
(512) H (234) or in dual twistors terminology to the point (13)w, one has 
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(AB (512) H (234)) (1345) 



(AB 12) (AB 23) (AB 34) (AB 45) (AB 51) 



AB 
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3.2 Chiral Integrals 

From the discussion of the pentagons, it is clear that there is a striking difference 
between a pentagon with a special numerator and plain scalar box integrals. Even 
though both kind of objects can be made pure integrals, each Schubert problem in 
the case of the pentagon has a single non-vanishing residue while in the boxes both 
solutions give rise to a residue. 

When an integral has the property that the residues associated to at least one 
of its Schubert problems are not the same, we say that the integral is chiral. The 
reason for the terminology comes from the fact that the two contours associated to a 
given Schubert problem are exchanged under parity (see section 2 for more details). 
This means that one can have chiral, pure, or chiral and pure integrals. 

At one-loop, one can have an even more especial class of integrals. When an 
integral has a numerator where at most one of the solutions to each Schubert problem 
gives a non-zero residue then we say that the integral is completely chiral. 

Let us give two more examples in this section. The first is the most general class 
of chiral pure pentagon integrals. This is an integral where only two of the five legs 
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needs to be massless. Moreover, it is clear that in order to write a special numerator 
the two massless legs cannot be adjacent. The claim is that the following family of 
integrals is (completely) chiral and pure. 
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In this case, the GRT can also be applied to show that all residues are the same. 
In order to show that the normalization gives unit leading singularities, identities of 
the form discussed at the end of this section are needed. 

Next, let us give a six-point two-loop example. Consider the following integral 



This integral has the structure of two of the general pentagon integrals joined by 
the all massive edge. Consider a residue of the full integral over C 8 which computes 
a residue of the pentagon on the left. The contour integral in Za and Zb is the 
same as before except that the normalization is different and therefore the residue 
is not equal to one. The residue must then be the ration of the two normalizations, 
i.e., (i j k I) j (i j CD). Plugging this in the integral over Z c and Z D we now find a 
properly normalized integral and therefore the remaining part of residue computation 
gives one. 

One might be tempted at this point to think that all completely chiral integrals 
are pure. In section 4, we describe in detail the example of a hexagon with a wavy 
line and a dashed line in the numerator. This integral is in fact completely chiral 
but it is not pure. 

3.3 Evaluation of Pure Integrals 

Evaluating integrals explicitly can be very hard and many techniques have been 
developed for this purpose. At one-loop, all integrals appearing in the standard 



./ 



k 




(AB i-l i) {AB ii+l) {AB j-l j) (AB j j+1) (AB CD) 
(CD (k-1 k k+l) f| (l-l I l+lj) 



{CD k-lk) {CD k k+l) {CD l-l I) {CD ll+l) 
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reduction techniques are known analytically. At higher loops, very few examples 
have been evaluated analytically. Many of our chiral pure integrals turn out to 
be completely IR finite and therefore their evaluation can be made directly four 
dimensions without any regulators. 

Consider the family of pentagon integrals discussed above. The evaluation of the 
integrals for generic j and k gives 
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{AB fl (j-ljj+l))(ijkk+l) 

{AB i-1 i) {AB ii+l){AB j-1 j) {AB j j +1){AB kk+l) ' 

= log (u^i-ij-i) log (w M _ MJ ) + Li 2 (1 - u jt fc.i-ij-i) + Li 2 (1 - u Ki _ hiJ ) 
- Li 2 (1 - Uj,fc,ij-i) - Li 2 (1 - Uij-i,k,i-i) + Li 2 (1 - 

(3.8) 

where 

= (ii+ljj+l)(kk+lll+l) , s 

Ui,j,k,i- (u + i jj+ i)(kk+lii+iy 1 J 

For special values of j and k the integral becomes IR divergent and a regulator is 
needed. We postpone this discussion to section 4. 

The reason for presenting the explicit form of the pentagon integrals is to note a 
general fact about pure integrals: The explicit evaluation of the integrals must be a 
linear combination of functions known as iterated integrals, such as polylogarithms, 
all with coefficient one. 

It is striking that the coefficients do not depend on kinematic invariants but 
this is a consequence of having unit leading singularities. This is the motivation for 
the terminology: pure integrals. Roughly speaking, the coefficients of the different 
polylogarithms are the leading singularities of the integrals. Having a pure integral 
ensures that no coefficient can depend on kinematical invariants. 

Once again, the hexagon with a wavy and a dashed line in the numerator given 
in section 4 will be an example of a completely chiral and IR finite integral which 
is not pure and its evaluation gives products of logarithms with different coefficients 
that depend on kinematic invariants. 

3.4 Example: One-Loop MHV Amplitude 

Up until know we have been studying integrals individually. This is a good point to 
actually use them to determine the full physical integral of the simplest set of ampli- 
tudes. These are one-loop MHV amplitudes. Historically, one-loop MHV amplitudes 
were the very first set of amplitudes to be computed for all n as a linear combination 
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of scalar box integrals [26]. It was found that the answer is very simple; an overall 
prefactor, proportional to the tree-level amplitude, and a sum over all one-mass and 
two-mass-easy box integrals with coefficient one, when properly normalized. In our 
modern terminology, the normalization was such that only pure integrals appear. 
It was realized that this form of the amplitude was not equivalent to the Feynman 
diagram amplitude as an expansion in the dimensional regularization parameter but 
it differs from it only at 0(e). In our language this is nothing but the fact that an 
expansion in terms of box integrals cannot possible reproduce the physical integrand 
of the theory as stressed a number of times already. 

Now that we have a set of chiral pure integrals, the natural question is how much 
more complicated the amplitude will look like if written in a form that matches the 
physical integrand. It turns out that the full integrand is stunningly simple 




(3.10) 



where the propagator (AB nl) is present in all terms. Note that not all integrals in 
the sum are chiral pure integrals. There are boundary terms which are box integrals. 
Consider for example j — % + 1. In this case the numerator cancels one of the 
propagators leaving us with the box. We give no derivation for this formula here and 
postpone a more detailed discussion to section 6. A final comment, even though the 
line (nl) seems especial, the amplitude is cyclic as it should be! 



4. Finite Integrals 

We have seen that the chiral integrals with unit leading singularities, naturally writ- 
ten in momentum-twistor space, provide a natural basis of objects to express the loop 
integrand. In this section we will see that they have another beautiful property- 
most such integrals are manifestly infrared finite. 

Let us begin by illustrating with a simple example. Consider a general 1-loop 
integral for 6 particles, which we can write as 

(ABX)(ABY) 



AB 



(AB 12) (AB 23) • • • (AB 61) 



(4.1) 
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Here X, Y are generic bitwistors. Of course, like almost all generic integrals with 
massless external legs, this integral is infrared divergent. Recall that the infrared 
divergences arise when the loop momentum / become collinear to a massless external 
momentum p a , i.e. when I • p a — > 0. The extra soft logarithmic divergence can be 
thought of as an even more special case of this situation, where the loop momentum 
becomes collinear to two consecutive momenta so that I -p a ,l -p a +i — > 0. In the dual 
co-ordinate space, the collinear divergence arises when the loop integration point x 
approaches one of the edges of the Wilson loop, connecting x a with x a+ i, and of 
course the extra soft divergence occurs when x gets close to both the lines (x a -i x a ) 
as well as (x a x a+ i), that is when it is close to the point x a itself. But again the IR 
divergence is fundamentally a collinear one, with the soft divergence being thought 
of as "double-collinear" . 

We can finally describe these IR divergent regions in momentum-twistor lan- 
guage. The collinear divergence associated with / • p a — > corresponds to the region 
where the line (AB) in momentum twistor space, associated to the loop integra- 
tion point, passes through Z a while lying the in the plane (Z a _iZ a Z a+ i). Note that 
this region is quite nicely parity invariant. Recall that in momentum-twistor vari- 
ables, parity is just the poincare duality, and exchanges the point Z^ with the plane 
W a i = (Z a ^iZ a Z a+ i)j naturally paired with Z a . Thus, the condition is that the line 
(AB) IJ passes through Z^, and also that the dual line (AB)jj passes through W a j. 

While a generic integral will indeed be IR-divergent, we see a simple way of 
getting completely IR finite integrals. If the bitwistors X, Y are chosen to have a 
zero in all the dangerous IR divergent configurations, then the integrals will be finite. 
This is very simple to achieve. For instance, let us choose X = (13) and Y = (46); 
we can write out the integral again as, 

6 1 

f (AB 13) (AB 46) (5 6 12) (2 3 4 5) 

J (AB 12) (AB 23) (AB 34) (AB 45) (AB56)(AB61) 1 ' ] 

AB 

4 3 

Let us check that the numerator has a zero in all the IR-divergent regions. Consider 
first collinearity with p 3 . We need to see what the numerator does when AB passes 
through Z 3 while lying in the plane (234). However, the numerator factor {AB 13) 
vanishes simply if AB passes through 1 or 3, regardless of whether or not it also 
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happens to lie in the plane (234). In this way, we can see that the collinear divergences 
with 1,3,4,6 are all killed by the numerator. Next, consider what happens when 
AB passes through 2, lying in the plane (123). Since AB lies in (123), it necessarily 
intersects the line (13), and therefore, {AB 13) = 0, regardless of whether or not 
(AB) also happens to pass through 2. A completely analogous argument holds for 
the collinear divergence associated with particle 5. 

Thus we see that with this numerator, all the regions with collinear divergences 
are killed by the numerator factors, and the integral is completely IR- finite! There 
are other choices for X, Y that will do the same job; our argument above also holds 
if one or both of the numerator factors (13), (46) were replaced by their parity- 
conjugates, (612) fl (234) and (345)0(561), respectively — changing one or more of 
the dashed-lines in (4.2) to wavy- lines. 

Now, these finite integrals are clearly chiral. And when the two numerators are of 
the same kind, they have, quite nicely and non-trivially, unit leading singularities. As 
usual, verifying by direct computation requires manipulating non-trivial sequences 
of 4-bracket Schouten identities, but the result follows much more transparently 
from an application of the global residue theorem to this integral. Consider for 
instance the GRT following from choosing j x = (AB 34), f 2 = (AB45), / 3 = (AB 56} 
and = (AB 61) (AB 12) (AB 23) . We have three different Schubert problems to 
consider, with the lines (34), (45), (56) combined with (61), (12), (23). Consider first 
the Schubert problem with the four lines (34), (45), (56), (61). This is a one-mass 
configuration, and it is easy to see that the numerator kills the solution where (AB) 
is the line (46), only leaving the solution passing through 5. Let us call this non- 
vanishing residue ^(34), (45), (56), (6i)- Similarly, for the Schubert problem with lines 
(34), (45), (56) and (12), the numerator kills the solution passing through 4 while 
leaving the one passing through 5; we can call this single non-vanishing residue 
?"(34), (45), (56), (12)- Finally, for the Schubert problem with lines (34), (45), (56), (23), we 
can see that both solutions — the line 35 as well the line passing through 4 — are killed 
by the numerator, so both of these residues vanish. The GRT then tells us that 

(0 + r (3 4),(45),(56),(61)) + (0 + r (3 4),(45),(56),(12)) + (0 + 0) = 

-> ?"(34),(45),(56),(61) = -?"(34),(45),(56),(12) (4.3) 

It is possible to repeat this argument for other GRT's, finding a sequence of 
2-term identities relating all the non-vanishing residues, showing that the integral 
is not only chiral but has unit leading singularities. Thus, we see in this instance 
something that can be checked also to be true for all other residues: the integral 
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is completely chiral; at most one of the two solutions to each Schubert problem are 
non-vanishing, and sometimes both vanish. 

Given that this integral has unit leading singularities, it is instructive to expand 
it in terms of boxes, which will then also have unit coefficients. This simple, finite 
momentum-twistor integral in fact expands into the sum of nine boxes: 




2 = { 



+ 



6 | 2 

5 4 

rj—K 

I 3 

5 6 
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The seemingly complicated combinations of a large number of boxes have been en- 
countered before in the computation of finite 1-loop objects, such as the NMHV ratio 
function [14,41-43] — the ratio function for the full superamplitude is simply defined 
to be 



n, 



1— loop 
n.k 



= A 



1— loop a tree a 1-loop 

n,k ^n,k ' • /l n,fc=2 ■ 



(4.4) 



Note that in the box expansion, every integral is individually IR-divergent, the 
IR-divergences only canceling in the sum. Moreover, the boxes themselves are not 
dual conformal invariant — again, only become dual conformal invariant in the sum. 
But since the hexagon in which we are interested is manifestly finite and dual confor- 
mal invariant 9 , we can evaluate it directly — for example, using Feynman parameter- 



9 In the literature on ratio functions, some authors have found what were claimed to be "finite" 
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ization directly without any regularization. A straightforward computation shows, 




2 = Li 2 (l - Ul ) + Li 2 (l - u 2 ) + Li 2 (l - u 3 ) + log(« 3 )log(ui) - y , (4.5) 



where the U{ are the familiar six-point cross-ratios 



Ml 



(12 34) (45 61) 
(12 45) (34 61) : 



«2 



(23 45) (56 12) 
(23 56) (45 12)' 



and 



«3 



(34 56) (61 23) 
(34 61) (56 23) 



(4.6) 



It is easy to find examples of integrals which are finite and chiral, but which 
do not have unit leading singularities. For example, changing one the 'dashed-line' 
numerator factor (AB 13) in the integral above to a 'wavy-line' [AB (612) f] (234)) 
will leave the integral finite and chiral, but spoil the equality of its leading singular- 
ities. Indeed, as it is also finite and dual-conformally invariant, the 'mixed' hexagon 
integral can also be evaluated without any regularization, and one finds that, 




2 = 



(AB (612)0 (234)) (AB 46) 



AB 



{AB 12) (AB 23) (AB 34) (AB 45) (AB 56) (AB 61) 



, (13^235) ) l0g(Ml) l0gM + ( (13f 5 )(5 4 613) ) l0gM l0g(Ml) + ( (1235) <lLl) ) l0gM l0gM ' 

In order for GRTs to yield the two-term identities necessary to guarantee that all 
the leading singularities are equal up-to a sign, the numerator must force vanishing 
residues for all but two Schubert problems. In the case of the 'mixed-numerator' 
hexagon integral, for example, GRTs can only be used to show that the coefficients 



combinations of boxes that did not end up being dual-conformal invariant. In every case, the 
combinations of boxes in question were not honestly IR-finite: the divergences from different regions 
of the integration contour canceling between each-other. Such a cancellation is is highly regulator- 
dependent, and is not very meaningful. 
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of the logarithms sum to zero: 



( (1234) \ / (6134) \ / (6123) \ 
V (1345) (1235); V( 13 45)(5613)y \ (1235) (3561) / 



0. 



(4.7) 



It is clear that these chiral momentum-twistor integrals with unit leading sin- 
gularities give us the simplest and most transparent way of talking about finite 
integrals. 

Just as a trivial example, the 6-point NMHV ratio function, which is typically 
written in terms of all 15 six-point box-integrals, with many /^-invariants as coeffi- 
cients, is given simply by 



integration 




} , (4-8) 



4 3 

x ([2 345 6] - [345 61] + [45 612]) 



where g : % \- > % + 1 acts on both the integrand and its coefficient. Also recall the 
definition of the R- invariants given in section 1, 



[ijklm] = 



£° |4 ({jklm)rjj + (klmi)r)j + (lmij)rj k + (mij k)rji + (ij kl)r) m ) 
(i j k I) (j k I m) (k I mi) (I mi j) (m i j k) 



(4.9) 
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5. One-Loop Integrands, Integrals, and Amplitudes 

As described in section 2.1, one can use elementary tensor-reduction to express any 
1-loop integrand in Af = 4 in terms of pentagon and box integrands. These of 
course would form a complete basis for any 1-loop integrand in Af = 4 SYM. How- 
ever, such a basis would necessarily include many integrands which are non-chiral 
(including all boxes), and which have non-uniform leading singularities; moreover, 
such a basis would allow for linear combinations of IR-divergent integrals to be ul- 
timately IR-finite and non-vanishing. But we saw in the last section that there are 
integrands — pentagons and hexagons with 'magic' numerators — which avoid all of 
these shortcomings, and these integrands closely mirror the leading singularities of 
MHV-amplitudes, suggesting that they may be well-suited to express amplitudes 
more generally. 

It is therefore natural to wonder if there exists a complete basis of 1-loop inte- 
grands involving only chiral, manifestly dual-conformally invariant integrands with 
unit leading-singularities, and for which no non-vanishing linear-combination of IR- 
divergent integrands is IR-finite. We'll see momentarily that the answer is affirma- 
tive, and extremely beautiful. 

Before trying to construct such a basis, however, we can gain some intuition 
about what to expect by assessing its size — that is, finding the dimension of the 
space of planar, 1-loop integrands. Recall that every n-point 1-loop planar integral 
can be written in the form 



When n = 5, the space of 1-loop integrands is just the space of bitwistors Y, which 
is six-dimensional — which explains how the complete 5-point 1-loop integrand could 
be constructed in [32] through the introduction of a single pentagon integrand. 

For n = 6, the most general integrand is a hexagon with {AB Yi)(AB Y 2 ) in the 
numerator. Now, each Yi is a 6-dimensional representation of SU4, and of course 
6 © 6 = 1 © 15 © 20. Ordinary multiplication being commutative, the antisymmetric 
part, the 15-component, clearly vanishes. By expanding each Yi into a basis of six 
simple bitwistors, it is easy to see that the trace component, 1, also vanishes, as 
(YiYj) = 0, when Y { is simple. Therefore, the space of 6-point 1-loop integrands is 
20-dimensional 10 . 

10 We thank Simon Caron-Huot for helpful discussions regarding this counting. 




(5.1) 



AB 
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More generally, it is not hard to see that the dimension, d, of 1-loop integrands is 
the same as the dimension of the space of symmetric n — 4-fold symmetric, traceless 
tensors of 6's of SU4, which is simply 

'=(i)+cr)- <-> 

Recall that box-integrands form a complete basis of parity-even integrands, and that 
there are precisely (j^j boxes, all of which are independent. Therefore, we may 

separate d in equation (5.2), according to d — d even + d dd with d even = (j^j an d 

d dd =(^ n 4 Once we have a basis of integrands which makes parity manifest, 
this will allow us to count the number of relations satisfied by (parity-odd) integrands. 

5.1 The Chiral Octagon: A Basis of One-Loop Integrands 

As we can see from equation (5.2), the number of independent integrands grows 
asymptotically like 0{n i ). In contrast, the number of chiral pentagons grows only like 
0(n 2 ). It is not hard to see that the simplest class of chiral integrands which number 
0(n 4 ) are the chiral octagons. As we will see presently, it turns out that chiral 
octagon-integrands indeed form an (over)-complete basis for all 1-loop integrands 
that satisfies all the desired criteria listed above. The most general chiral octagon 
integral is given by, 



I 8 (i,j,k,l) = If 3Y~ J for i<j<k<l<i (5.3) 




k 



(ABij)(AB{j-lj j+l)f](k-l k k+l))(ABkl)(AB(l-l I Z+l) PI i 
{AB i-1 i) {AB ii+l) (AB j-1 j) (AB j j +1) (AB k-1 k) (AB kk+l) {AB l-l I) (ABll+1)' 

AB 

Notice that parity acts according to P : Ig(i,j,k,l) 1— >■ I$(j,k,l,i), making it trivial 
to define parity-even/parity-odd sectors: 



C en/ ° dd (^i>M) = h(hj,k,l)±I 8 (j,k,l,i) 
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Figure 3: The possible degenerations of the general octagon integrand. 

i 




k 



Notice that these octagon integrands are well-defined for any distinct set of 
indices {i,j,k,l}, including those for which the 'octagon' degenerates into lower- 
polygons. For example, when I = k + 1, the extra (duplicated) propagator in 
equation (5.3), (AB k k+1) is cancelled by the dashed-line term (AB k I) — > (AB k k+1) 
in the numerator. A complete sampling of degenerate 'octagons' is illustrated in 
Figure 3. One important advantage of this presentation is that all but the most de- 
generate of the octagons are manifestly finite. Indeed, only the pentagons — octagons 
of the from I s (i, i+2, i+3, i+4) — and the lower hexagons in Figure 3 — octagons 
of the form I 8 (i,i+l,i+2, j) — are IR-divergent. Specifically, we have the following 
separation into manifestly IR-finite and IR-divergent basis integrands. 




It is worth emphasizing that the only IR-finite combinations of IR-divergent integrals 
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in this basis are parity-odd, which automatically vanish upon integration. Further- 
more, as discussed above, because the criterion for local divergences in the region of 
integration is itself parity-invariant, 11 parity-odd combinations of integrands are in 
fact manifestly locally finite. 

Parity-symmetrizing, and parity anti-symmetrizing, it is clear that there 2 ^ ™ 

octagon integrands, evenly split between parity-odd and parity-even. As we described 

n — 1 



above, among the parity-odd combinations of integrands only y ^ j are linearly- 
independent, so the octagon basis is strictly over- complete, but there are only non- 
trivial relations among integrands in the parity-odd sector. 

5.2 Integration of Manifestly-Finite Octagons 

It is not hard to directly evaluate the general octagon integral h(hj, k, I)- Consider 
for example the case J 8 (3, 6, 9, 12) for which all indices are separated by at least 3, 



(5.6) 




9 



Because of the numerator factors, the only non-vanishing leading singularities of this 
integral involve cutting at most one of each the pairs of lines {(23), (34)}, {(56), (67)}, 
{(89), (9 10)}, and {(11 12), (12 1)}. Therefore, this integral's box-expansion is sim- 
ply the (manifestly-finite) sum of 16 four-mass box integrals. One disadvantage 
with this presentation of the integral, however, is that the four-mass box integral 
logarithmically-diverges when any of its four massive corners becomes massless, and 
yet we saw above that the general octagon remains manifestly finite upon many such 
degenerations. 

Because of this, we are motivated to replace the four-mass box function with a 
new function that is free of any divergences over the physical domain of cross ratios. 
Letting A4(u,v) denote the familiar four-mass box integral 12 — a symmetric function 



11 Recall that a local, IR-divergence develops in the region of integration when the line (AB) 
passes through a point Zi while simultaneously lying on the plane Zi Zi + \). 

12 If we denote the massive, incoming four- momenta of the box by K\ , K 2 , if 3 , and K4, and define 
the canonical Mandclstam variables s = {K\ + K 2 ) 2 and t = (K 2 + K3) 2 , then we are using u and 
v to denote the cross ratios K 2 K 2 /(st), and K\K\j{st), respectively. 
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in the two cross-ratios — then let us define the following 'modified four- mass' function 



A 4 (i,j,k,l) = Ai{ui^k,h^j,k,hi) ~ 2^°&( Ui d,k,i)\°g( u j,k,i,i) 



(5.7) 



where 



and 



A 4 (u, v) = Li 2 (l - a + ) - Li 2 (l - ot_) + - log(f) log(a + /ot_), 



2u 



ot± 



(5.8) 



(5.9) 



1 + u — v + — u — v) 2 — Auv 

here, we have used the four indices {i,j, k, 1} to signify the (generally time-like 
separated) spacetime points corresponding to the lines (j (kk+1), and 

(ll+l) in twistor space, which together define the cross-ratios 

__ (ii+ljj+l)(kk+lll+l) _ {j j+1 kk+1) (ll+l 

Ui,j, k ,i= (u + i jj+ i)(kk+iii+i) u ^,i,i = (ii + ikk + i)(ii + ijj+iy 

(5.10) 

The principle distinction between A 4 (z',j, k,l) and the more familiar four-mass 
box function is that A 4 (i,j, k,l) remains finite even when many of the spacetime 
points become null-separated (or even become identified). In particular, 



lim ( Ai(i,j,k,l) J = Li 2 (l -Ujk,n) and lim ( A 4 (i, j, k, I) ) = Li 2 (l -u ijk! i)- 

u i,j,k,i->0 V / Uj,k,i,i->0 V / 

(5.11) 

Of course, if we use A 4 's to represent Is(3, 6, 9, 12), for example, then each four- 
mass box will contribute a 'log-log'-term. It may be worried that this will greatly 
clutter the final expression, but this turns out to not be the case: taken together, 
these 16 additional 'log- log' terms combine into a single such term. 

With this new function, the general octagon integral — together with all its 
degenerations — becomes extremely simple. Explicitly, the general octagon I s (i, j, k, I) 
integral is given by, 




log (Uik-lki-l) \0g(Uji-iij-i) 

+A 4 (i,j,M) -A 4 (i,j,M-l) -A 4 (i,j,fc-l,Z) +A 4 (i,j,fc-1,/-1) 
-Aili,j-l,k,l) +A 4 (z,j-1,M-1) +A 4 (i,j-l,k~l,l) -A 4 (z,j-l,fc-l,/-l) 
-A 4 (i-l,j,k,l) +A 4 (i-l,j,k,l-l) +A 4 (i-l,j,k-l,l) -A 4 (i-l,j,k-l,l-l) 
[+A 4 (z-1, j-1, k, 0-A 4 (z-l, j-1, k, Z-1)-A 4 (z-1, j-1, k-1, 0+A 4 (z-l, j-1, k-1, l-l) 
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k + l 



k + l 



Although admittedly lengthy, this expression can be considerably compressed in a 
way which helps illustrate the relative signs appearing in the formula above, 

h(i,j,k,l) =log(u ifc _ lfci _ 1 )log(u J -,_ 1Jj _ 1 )+ (-l) {ai+a2+a3+a4) ^S-vi,J-V2,k-a 3 ,l-a 4 ). 

^£{1,0} 

(5.12) 

We can see how the modified four-mass function A 4 (i,j, fc,Z) helps to make all 
of the octagon's degenerations manifest by looking at a few examples explicitly. For 
example, consider the 8-point octagon 7 8 (2,4, 6,8); in this case, only 20 of the 34 
cross ratios which play a role in the general answer are non-vanishing, converting 
virtually all the generalized four-mass functions A 4 's into Li 2 's. 

2 

3 

log(W2,5,6,l) lo g( M 4,7,8,3) 
+ A 4 (2,4, 6,8) -Li 2 (1 - M 4i6i 7 i2 )-Li2 (1 - M 2 ,4,5,8)+Li 2 (1 - M2,4,5,7) 
-Li 2 (1 - M 8)2 , 3 ,6)+Li 2 (1 - M 3 ,6,7,2)+Li 2 (1 - M 8 ,2,3,5)-Li 2 (1 - U 7 ,2,3,5) 
-Li 2 (1 - M 6 , 8 ,l,4)+Li 2 (1 - M 4 ,6,7,l)+Li 2 (1 - Wl,4,5,8)-Li2_(l - ^1,4,5,7) 

+Li 2 (1 - 1*6,8,1,3) -L12 (1 - W3,6,7,i)- Li 2 (1 - "5,8,1,3)+ A 4 (l, 3, 5, 7) 

Even more simplification occurs for the degenerate 'octagons.' Consider for 
example the general finite heptagon integral, given by fc, fc + 1), 

log (Ui.fc-i.fc.i-i) log (Uj,fc,fc+lj-l) 

+Li 2 (1 - u jAk+hi ) -0 - A 4 (i, j, fc-1, fc+1) +Li 2 (1 - u jtk _ 1Ai ) 
-Li 2 (1 - Uj-i,fc,fc+i,i) +0 +A 4 (i, j-1, fc-1, fc+1) -Li 2 (1 - Uj-i >k -i,k,i) 
-Li 2 (1 - Uj,fc,fc+i,i_i) +0 +A 4 (i-1, j, fc-1, fc+1) -Li 2 (1 - u j>k -i,k,i-i) 
+Li 2 (1 - Uj-i^.fc+i.j-i) -0 -A 4 (i-1, j-1, fc-1, fc+l)+Li 2 (1 - w^x^x^x) 

Here, because A 4 (i, j, fc, fc) = Li 2 (0) = 0, four of the contributions vanish, and eight 
of the modified four-mass box functions simplify to simple Li 2 's. 

The final class of finite, degenerate octagons are the hexagon integrals — octagons 
of the form i + 1, fc, fc + 1), 






i + l 



log (u ijk -l,k,i-l) log («i+l,fc,fc+l,i) 

+Li 2 (1 — u i+ i ikjk+lii ) — — Li 2 (1 — Mfc+i,i,i+i,fc-i) +Li 2 (1 — Uj+i,fc-i,fc,i) 

-0 +0 +0 -0 

-Li 2 (1 - Ui+i,fc,fc + i,i_i) +0 +A 4 (i-1, i+l, fc-1, fc+l)-Li 2 (1 - Ui+i,fc-i,fc,i-i) 

+Li 2 (1 — Wi,jt,fe+i,j-i) —0 — Li 2 (1 — , Ufc+i,i-i,j,jfc-i) +Li 2 (1 — Uj^-i^^-i) 
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Just as for the case of the 8-point octagon integral, the general hexagon integral 
simplifies considerably when potentially-massive corners become massless. As a final 
illustration, let us see how the general formula for the octagon given above directly 
yields the result quoted in section 4 for the 6-point hexagon integral which played 
such an important role in the 6-point NMHV ratio function: 



log (113,5,6,2) log (114,6,1,3) 
+Li 2 (1 - ^4,6,1,3) -0 -Li 2 (l) +Li 2 (l) 
-0 +0 +0 -0 

-Li 2 (l) +0 +Li 2 (1 - M 2 , 4 ,5,i)-Li 2 (l) 

+Li 2 (l) -0 -Li 2 (l) +Li 2 (1 - u 3i5)6)2 ) 




= Li 2 (l - 114,6,1,3) + Li2(l - ^2,4,5,1) + Li 2 (l - 1*3,5,6,2) 
3 + log(u 3) 5 >6)2 ) log(u 4) 6,i,3) - 2Li 2 (l). 

5.3 Application: the NMHV One-Loop Ratio Function 

As should be clear from the previous subsection, any IR-finite object such as the 
ratio function will be manifestly finite when expanded in the basis of octagon in- 
tegrands. Moreover, since the formula for the completely general octagon integral, 
equation (5.12), is free of discontinuities for all the IR-finite degenerations of the 
octagon, any finite 1-loop integrand expressed in the basis of octagons directly trans- 
lates into a function that is manifestly dual-conformally invariant. 

A very important, manifestly finite function associated with 1-loop scattering 
amplitudes is the ratio function, 



1 — loop 



aI— loop jtrcc jl— loop 

^n,k ^n,k ' • /i ra,fc=2 ■ 



(5.13) 



The most trivial example must be the 5-point 1-loop NMHV ratio function. Expand- 
ing into the basis of octagons, the integrand is easily seen to be given by 

/ 2 2 \ 



^5,3,1 = [1,2,3,4,5] 




+ cyclic. (5.14) 



v - - - - / 

Being a parity-odd combination of pentagons, the ratio function is locally free of any 
divergences at the level of the integrand and is therefore manifestly finite — of course, 
being parity odd, it also vanishes upon integration. 
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A less trivial example, and one which we quoted in section 4, is the 6-point 
NMHV 1-loop ratio function. In section 4 only the parity-even contribution to the 
ratio function was described; the full integrand is given by, 



7^6,3,1 




+ 



5 4 

x \ ([1,2, 3, 4, 5] + [1,2, 3, 5, 6] + [1,2, 3, 6, 4]) 



/ 1 



+ 




3 — 6 




\ 



+ 



/ 1 




2 3 



V 5 4 5 4 / 

x| ([1,2, 3, 4, 5] -[1,3, 4, 5, 6]) 

+ cyclic. 



\ 



5 4 5 4 

xMl,2,3,4,6] 



2 3 




V 6 5 6 5 / 

x|([l,2, 4,5,6] + [1,3,4,5,6]) 

(5.15) 



Of course, only the first term in equation (5.15) is non-vanishing when integrated 
along a parity-invariant contour, reproducing the formula given in equation (4.8). 

The general formula for the n-point NMHV 1-loop ratio function integrand nicely 
separates into a part which is parity-odd, and another which involves only manifestly 
finite integrands. In order to best capture the ratio function succinctly, let us intro- 
duce one small bit of notation and define 

O-iJ} {l,k} 

+ {k,. ..,/>]= Yl E M>*1; ( 5 - 16 ) 

J={i+l,i+2} K={k,k+1} 



for example, 



[1, {2, 3, 4}, {5, 6, 7}] = [1, 2, 3, 5, 6] + [1, 2, 3, 6, 7] + [1, 2, 3, 7, 5] 

+ [1,3, 4, 5, 6] + [1,3, 4, 6, 7] + [1,3, 4, 7, 5]. 



(5.17) 



Notice that the two-index J ranges over all consecutive pairs between i + 1 and j 
inclusively, while the two-index K also includes a non-consecutively-ordered 'wrap- 
ping' term. With this notation, it is very easy to write the n-point NMHV ratio 
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function integrand: 

KT P = \ E ih{i+h...jh{k,...,m(ij,k,i) 

i<J<k<l<l (5.18) 

i<j<k<l<m<i 

Notice that while the first term in equation (5.18) appears to include divergent 
'octagons', only the finite octagons have non- vanishing coefficients. For five-particles, 
for example, the coefficient of the octagon 7 8 (1,3, 4,5) from equation (5.18) would 
be [1, {2, 3}, {4, 5}] = [1,2, 3, 4, 5] + [1, 2, 3, 5, 4] = 0. 

Combining formula (5.18) with the analytic form of the general octagon inte- 
gral given in equation (5.12) immediately yields a concise, analytic, manifestly dual- 
conformally invariant, and manifestly-cyclic form of the 1-loop ratio function for any 
n. 

Let us close this section by given another explicit example. The 7-point NMHV 
1-loop ratio function is straightforwardly found to be, 

U 1 ^ = [1, {2, 3}, {4, 5, 6}]/ 8 (l, 3, 4, 6) + [1, {2, 3}, {4, 5, 6, 7}]/ 8 (l, 3, 4, 7) + cyclic, (5.19) 

where 

J 8 (l, 3, 4, 6) = | Li 2 (1 - "1,3,4,5) +Li 2 (1 - "2,4,5,1) + Li 2 (1 - "4,6,7,2) +Li 2 (1 - "7,2,3,5) 
- Li 2 (1 - "2,4,5,7) -Li 2 (1 - "3,6,7,2) -Li 2 (1 - "4,6,7,3) -Li 2 (1 - "6,1,2,4) 

+ l0g ("1,3,4,7) log ("3,5,6,2) }; 

I 8 (l, 3, 4, 7) = I log ("1,3,4,7) log ("3,6,7,2) - Li 2 (l) -Li 2 (1 - 1*1,3,4,5) -Li 2 (1 - "4,6,7,2) 

+ Li 2 (1 - 1*1,3,4,7) +Li 2 (1 - 1*3,6,7,2) +Li 2 (1 - "4,6,7,3) +Li 2 (1 - "6,1,2,4) }• 

(Here, we have not neglected an overall factor of |: like in the case of the 6-point 
ratio function — the summand in equation (5.18) includes exactly two copies of each 
term; but this is not generally the case for higher-n). 
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6. Multiloop Amplitudes 



In this section, we introduce a new strategy for finding local representations of loop 
integrands. The idea is closely related to the leading singularity method, but the 
philosophy differs in some important ways. In particular we will not be guided by sys- 
tematically trying to match all the leading singularities of the integrand. Instead, we 
will look at a simple subset of leading singularities defined for generic, large enough 
number of particles — no "composite" leading singularities will be considered. We 
will then find a natural set of pure integrals designed to match this subset of leading 
singularities. We will find that boldly summing over all such objects miraculously 
suffices to match the full integrand! In particular, while the pure integrals are mo- 
tivated for a large-enough generic number of external particles, their degenerations 
nicely produce all the needed lower-point objects as well. 

This method is heuristic — we do not yet have a deep understanding for why the 
miracles happen. However we have used this strategy successfully to find stunningly 
simple expressions for the integrands of all 2- and 3-loop MHV amplitudes as well 
as all 2-loop NMHV amplitudes, and have checked that the results are correct by 
comparing with the form obtained from the all-loop BCFW recursion. 

We will begin by illustrating this strategy by going back to 1-loop integrands, 
which will motivate structures for 1-loop integrands different from the ones we en- 
countered in section 3. For the MHV integrand, this new form coincides with one 
of "polytope representations" discussed in [9]. We will then use this discussion as a 
springboard to our treatment of 2- and 3-loop integrands. 

6.1 New form for the MHV 1-loop Integrand 

Let's begin by going back to the MHV 1-loop integrand, and motivate a new form 
for it inspired by straightforwardly matching its leading singularities, associated with 
the familiar two-mass-easy colored diagrams 



i . . i 




corresponding to cutting the propagators 

(AB i-1 i) (AB i (AB j-1 j) (AB j 



(6.1) 



The amplitude has unit leading singularity for the first solution of the Schubert 
problem (AB) = (ij), and vanishing leading singularity for the second solution 
whereAB = (i-lii+1) f)( j-1 j j+1). We would like to build the integrand out 
of objects that have exactly this property. To beat a dead horse yet again - 
it is obvious that the two-mass-easy box does not do this job because it is not 
chiral. The easiest way to do this is to simply insert a factor in the numerator, 
(AB (i-lii+1) H O'-ljJ+l)), that kills the "wrong" leading singularity. For correct 
little-group weights, we add a factor (ABX) in the denominator, where X is an 
arbitrary bitwistor, and look at an object of the form 



which is just the pentagon already familiar from section 2, where the local propagator 
(ABnl) has been replaced by (ABX). We denote this graphically as 



Note that there is in general no significance to the presence of the legs adjacent to X 
in this picture. We draw it in this way because in the special case where X — (k k+1), 
the legs adjacent to X are identified with k,k + l. 

Now consider the Schubert problems associated with cutting four physical propa- 
gators. By construction this object has vanishing leading singularities on the "wrong" 
solution, and can easily be seen to have unit leading singularity on the "right" one. 
Summing over all the indices % < j — with \i — j\ > 2 corresponding to the two-mass 
easy colored graphs — produces an object matching all the physical leading singulari- 
ties of the amplitude. Naively this should give us the integrand, but there is a catch: 
each term also has "spurious cuts" where (ABX) is one on the cut propogators. 
Indeed, the sum we just described does not match the integrand. 

However some wonderful magic happens: the sum over all indices % < j, including 
a "boundary term "with j = i + 1, which is not included in the sum over colored 
graphs, does reproduce the amplitude! We have 



(AB (i-lii+l)[\(j-ljj+l))(Xij) 



(6.2) 



(AB i-1 i) (AB i i+1) (AB j-1 j) (AB j j+1) (AB X) 



J 




(6.3) 



A 



1— loop _ 
MHV ~~ 



E 



(AB n (j-ljj+l))(Xij) 



(6.4) 



i<j 



(AB X)(AB i-li)(ABii+l)(AB j-1 j)(AB jj+1) 
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or in the form of picture 



a1— loop 

•^MHV ~~ 




(6.5) 



This form is manifestly cyclic but has spurious (ABX) poles term-by-term. The sum 
is however independent of X. If we choose X to correspond to one of the external 
point X = (kk + 1), all the poles are manifestly physical but the formula is not 
manifestly cyclic invariant. 

As mentioned above, this expression follows from a simple "polytope" interpreta- 
tion [9]. The local formula given in [1] is obtained by choosing X = kk+1, summing 
over all k and dividing by 1/n. The similar expression in [5] corresponds to setting 
X = Jqo where 1^ is infinity twistor. 

Let us look at the "boundary term" where j — % + 1 in more detail-using 
ii+l i+2) {AB ii+l) = {AB fl (ii+li+2)), we can see that it is just 

a (spurious) box 



It is instructive to explicitly understand the purpose of this boundary term in this 
simple example, since the same phenomenon will occur in all the rest of our examples 
in this section. Let us return to our most naive ansatz, summing only over the 
pentagons associated with the colored graphs. Each of the spurious cuts involving 
{ABX), such as 



is shared by two pentagons e.g. hj-i and Tjj. For generic terms in the sum, these 
cuts cancel against each other in pairs. However, in the limiting cases when j = i+2 
(or j = i-2) the quad-cut is shared by 1^+2 and Jj-i^+i but there is no cancelation 
between them because the non-vanishing leading singularities occur for two different 
solutions of the Schubert problems. The spurious box of (6.6) precisely has non- 
vanishing leading singularities for these two Schubert problems and completes the 
cancelation of all {ABX) poles, ensuring the full sum is independent of X. It is quite 
remarkable that the "new" object needed to fix the leading singularities and match 
the amplitude is simply a degeneration of the pentagon itself. 



{i-1 ii+l i+2) {X ii+l) 



(6.6) 



{AB X) {AB i-1 i) {AB ii+l){ABi+l i+2) 



{ABX){AB i-li){AB ii+l) {AB j-1 j) 



(6.7) 
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In our remaining examples, we will not delve into understanding the details of 
how all leading singularities match. We will instead take a class of leading singu- 
larities as a guide for the local integrals to consider, and sum over all the relevant 
objects, including boundary terms that do not directly correspond to any of the lead- 
ing singularity pictures that motivated the construction of the objects to begin with. 
These formulae are then verified by comparing with the integrand as computed by 
BCFW recursion. 

Let us finally note a very pretty property of eqn.(6.4): for generic X, all the 
pentagons in the double sum are manifestly manifestly IR finite. This ceases to be 
true if we make the special choice like X = (12), since the diagrams with i = 2 or 
j = n have an additional massless corner which is not controlled by the numerator. 




6.2 NMHV 1-loop integrand revisited 

We proceed to use the same strategy to determine a local expression for the NMHV 
1-loop integrand, which will yield a quite different form than we obtained in section 
3. We again start with the colored graphs for leading singularities. There are two of 
them for NMHV amplitudes: 




Unlike the MHV case where the non- vanishing leading singularities were "1", 
here the non- vanishing leading singularities are the i?-invariants. The goal is to find 
objects with non-vanishing support on the same Schubert problems as the amplitude, 
and decorate these with the appropriate /^-invariants to get a nice ansatz for the 
integrand. 
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The first colored graph correspond to 2-mass easy Schubert problems and have 
the same structure as the MHV case. The leading singularity is just the tree-level 
amplitude appearing in the upper-left corner of the figure, A^^ K y(j, j+1, • • • ^-1, «)■ 
Thus we expect to have objects in the integrand of the form 




(6.10) 



Finding an object matching the physical leading singularities of the second class 
of colored diagrams is a more interesting exercise. The cut propogators are 

(AB i-l i) (ABi i+1) (ABj j+l)(AB k k+1) (6.11) 

The leading singularities vanish for the solution AB = (i-l H (j J +1) (i— 1 i i+1) f] (k k+1), 
while for AB = (ij j+1) f](ik k+1) the leading singularity is [i,j,j+l, k, k+1]. 
Let us consider objects of the form 

i 

[ N(i,j,k) 

J (AB X) (AB i-l i) (AB ii+l)(AB j j+1) (AB k k+1) 

AB 

We are searching for a numerator supported on the same leading singularities as 
the amplitude. In addition it should also have unit leading singularity on all other 
spurious quad-cuts. The reason is that the spurious cuts must cancel in a sum over 
terms; since the integrals are multiplied by different /^-invariants, the only way this 
can happen is through residue theorem 6-term identities between the i?-invariants. 
For instance the spurious quad-cut 

(ABX) (ABi i+1) (AB j j+1) (ABk k+1) (6. 12) 

is shared by six different integrals J iy - )fc , h+i-j^, Ij-i,k, Ij+i-,i,k, h-^j and 4+i ;iJ that 
are multiplied by six different residues. There is a 6-term identity relating them 

[i,j,j+l,k,k+l] + [i+1, j, j+1, k, k+1] + [j,i,i+l,k,k+l] 

+ [j+l,i,i+l,k,k+l] + [k, i, i+1, j, j+1] + [k+1, i, i+1, j, j+1] = 
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which can only possibly be of help in canceling spurious cuts if the integrands they 
multiply have support on the same Schubert problems, with unit leading singularities. 

There is one final guiding principle for determining the structure of the numerator 
N(i,j, k). The topologies occurring in (6.10) are the same as for the MHV amplitude, 
while the second class of integrals is "purely" NMHV-like. Since IR-divergences are 
universal, it would be nice if the IR-divergences could be completely isolated in the 
MHV-like topology. We should then try to choose the numerator N(i,j,k) to be 
strictly finite. It would be nice if these integrals could be chosen to be manifestly 
finite. The only divergence in (6.12) can come from the Zj-corner, i.e. the region when 
AB crosses point Zi and lies in the plane (i—\ In order to control this region 

the numerator should be of the form N = (AB fl (• • • ))• Combined with 

the unit leading singularity constraint, the form of the numerator is fixed completely: 

N(i,j,k) = (AB(i-lii+l)r\^ij, k ) (6.13) 
with £j,j,fc a special plane defined according to 

E iJtk = l -[{jj+l (ikk+l)f\X) - {kk+l (ijj+l)r[X)] (6.14) 

This is in fact the only choice we could have made consistent with little group weights 
and the desire to treat the j, k indices symmetrically. We will denote this by, 



i 




(ab (a»+i)ns i# ) 

(AB X) (AB i-1 i) (AB ii+l)(AB j j+l)(AB k k+1) 



With these objects in hand, we once again brazenly sum over all ranges of indices, 
including "boundary" terms with j = i±l not directly associated with colored graphs 
for leading singularities. The same magic happens as we saw in the MHV case — this 
sum agrees with the 1-loop NMHV amplitude as computed by BCFW recursion, and 
we find, 
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i 




(6.15) 



Note also that as in the MHV case, the only IR divergent integrals are in the 
boundary terms. The (generically) finite integrals for ly^ are given by 

Ii,j,k = -Li 2 (1 - «i) - Li 2 (1 - u 2 ) + Li 2 (1 - u 3 ) + log (w 4 ) log (u e ) + A 4 (u 4 , v 4 ) + A 4 (u 5 , v 5 ) 

where the cross ratios are defined as: 

(ii+ljj+l)(Xi-li) _ (ii+lX)(kk+li-li) _ (ii+lj j +1) (k k+1 i-1 i) 

Ul = (ii+lX)(j j+li-li) ' U2 = (ii+lkk+l)(Xi-li) 7 ^ 3 = (i i+1 k k+l) (j j+1 i-1 i) ' 
(Xkk+l){i~lijj+l) _ (jj+lX){kk+li-li) _ (j j+lX){kk+lii+l) _ (X k k+l) {ii+lj 
(Xi-1 i)(kk+lj j+1) 1 Vi = (jj+lkk+^iXi-li) 1 U5 = (jj+lkk+l)(Xii+l)' V5 = (Xii+l){kk+l jj+l) 

Finally, let us examine the 1-loop NMHV ratio function 

-r,l-loop _ j 1-loop _ 4 1-loop 4 tree (c -i p\ 

NMHV ^NMHV "^MHV ' "^NMHV l u - iu i 

Comparing the expressions 6.4 and 6.15 we can see that the ratio function 
has the same form as NMHV amplitude, except that in the first sum we have 

^NMHv(*> • • • J- 1 ' j) -^NMHV in stead Of just ^nmhv(^ ■ ■ ■ The man " 

ifest finiteness is obvious. The only divergent integrals are in the boundary term j = 
i — but their coefficient is given by A^ H y(i,i+l, . . — ^nmhv(1> • • • > n ) = 

-^nmhv(*i i+1, • • • i-1) — ^•nmhv(-'-) • • • > n ) = ®' Therefore, the ratio function can 
be written only using manifestly finite integrals. 

6.3 The 2-loop MHV amplitude and its logarithm 

Now we turn to the 2-loop case. First we reproduce the MHV amplitude presented 
already in [1] and in addition we will write an expression for the log of the amplitude 
given in an interesting form in terms of non-planar diagrams. 
We again start with colored graphs, 
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There are more types of graphs in comparison to 1-loop where we had only 
boxes. In addition to two glued boxes (also referred to as "kissing boxes") we have 
other topologies — pentaboxes and double-boxes. They represent cutting the internal 
(ABCD) propagator once and twice respectively, the latter case corresponding to 
"composite" leading singularities. 

Let us concentrate on the first graph. It looks like a "squaring" of the 1-loop 
cuts with appropriate ranges for indices. And in fact, the AB part and CD part of 
the integral are independent, ie. in order to realize the octa-cut of the first colored 
graph, we need to set AB = ij or AB = (i-lii+l)f] (j— and CD = kl or 
CD = (k—1 k k+1) Pi (i—lii+l). Together we have four possible combinations. The 
amplitude (as we see from the colored graph) has support just on one of them (AB) = 
ij and CD = ki while for all other it vanishes. It means that the numerator must 
vanish whenever AB = (i-1 i i+1) f| (j-1 j j+l) or CD = (k-lk k+l)f](£-l££+l). 
This motivates us to start with an integral of the form 



Jk. 

X 



(ab j/ l/z-nni./ !././• D);.v/./) 

(AB X) (AB (AB ii+l)(AB j-1 j) (AB j j +1) 
(CD (k-1 k k+1) (1 (l-l l l+l)){Y k I) 

(CD Y) (CD k-lk) (CD k k+1) (CD l-l I) (CD I l+l) 



x 



which has exactly this property. However, there is a better candidate. Instead of 
adding (ABX) and (CDY) in the denominator, we can add directly the internal 
propagator (ABCD). That allows us to write two numerator factors exactly as we 
need. Therefore, we consider, 
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(AB (/ l//-nni,/ !././• U)(/./ ^-/) 
(AB i-1) {AB i i+1) (AB j-l j) (AB j j+1) (AB CD) 

(CD (k-1 k k+l) fl (l-l I l+l)) 
X (CD k-lk) (CD k k+l) (CD l-l I) (CD I l+l) 



Of course, this integral has also many other cuts - both composite and non-composite 
- that involve the propagator (ABCD) , and we have to match other colored graphs 
in (6.17) as well. However, just as in our 1-loop examples, simply summing over all 
indices with a planar ordering reproduces the full amplitude as a cyclic sum over just 
one integral topology: 



a 2— loop 
•^MHV 



i<j <k<l<i 




(6.18) 



The "boundary terms" in this case occur for for j = i+1 and/or I — k+l. In these 
cases the numerator exactly cancels one of the propagators, leaving us with 13 





(6.19) 



Log of the amplitude 

Finally, we give an interesting new expression for the logarithm of the amplitude, 
using a non-planar sum of the same set of objects. At 2-loops, the log of the amplitude 
is 



[iog-4r ioop 



jmhv 



a 2— loop 
•^MHV 



1 (a 1-loop V 

2 y MHV J 



(6.20) 



13 This simplification was missed in [1], and the 2-loop MHV integrand was presented as a sum 
over three terms. We would like to thank Johannes Henn for pointing the simplification out to us. 
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A beautiful expression for the log of the amplitude is made possible by the 
existence of a simple relation between the sum of 1-loop square and 2-loop diagrams: 




(6.21) 



The left-hand side is just (Al A ^ > y P ) 2 while the right-hand side contains not only 
the planar diagrams present in Amhv P but also non-planar graphs when for example 
i < k < j < I. In fact, all planar graphs are equal to 2.4^ 1E ^ P while all non-planar 
graphs give us the log of the amplitude in the form 

j k 




[log^MHV P = -^ X \ f X (6-22) 

i<k<j <l<i 

i I 

The formula found in [27] is the 4pt version of this expression. 

Note that naively, all these integrals are IR finite because each individual 1- 
loop sub-integral is just a finite pentagon(which can not shrink to a box due to the 
restriction j ^ i + 1 and I ^ k + 1). However, the criteria for finiteness we described 
in section 4 applies to planar integrals, while the log contains non-planar terms which 
can be IR divergent. 

Let us focus on the piece of the integrand of the form 

(ABX) 1 (CDY) 

(ABi-li)(ABii+l) (ABCD) (CDj-1j)(CDjj+1) 1 ' ' 

Here X controls the IR divergence of the region where the line AB intersects point 
Zi and lies in the plane (Zi-iZiZt^), just as Y does for CD sector. However, 
if % — j then AB and CD intersect in the point i and the propagator (ABCD) 
vanishes. Therefore, finiteness of the 1-loop sub-integrals is not enough. We need an 
extra condition that regulates this joint divergence. It is not hard to see that unless 
(XY) = 0, a (mild) IR divergence remains. 
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As a result, we can find that almost all integrals in (6.22) are finite except for 
the class of diagrams: 

i + 1 i 




(6.24) 



In this case X = (i-2i-li) C\(ii+li+2) and Y = f| (J-l j j+1), so 

(XY) 7^ 0. However the divergence is mild, as observed in the 4-pt result of [27]. 

6.4 2-loop NMHV amplitude 

We move on to present the integrand for all 2-loop NMHV amplitudes. The 6- and 
7- point integrands were presented in [1], by expanding the BCFW result into a basis 
of pure integrals. The parity-even part of the 6- point integrand was presented using 
standard (dual) space-time variables in [8]. Here, instead of a brute-force expansion 
into a basis of integrals, we follow the same strategy outlined above, obtaining results 
vastly simpler than those presented to date, which also generalize to all ra.Now 

Let us first start by drawing the colored-graphs that contribute for general 2-loop 
NMHV amplitude that do not cut the internal propagator (ABCD). 



j k j k j k 




[ijj+lkk+1] [ijj+lkk+1] [ijj+lkk+1] 
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Below each colored graph, we have indicated the leading singularity below each. 
Notice that the coefficient ^nmhvC?' ■ k;l, .. .i) is the same function as an ordinary 
tree amplitude with particles labelled (J, . . . ,k;l, . . . i) where k, I and i,j are both 
treated as if they were adjacently-labelled. 

The idea is again to find a set of integrals that each individually have the same 
leading singularities as the amplitude on a given set of octa-cuts. The first step 
is to realize that the octa-cuts on the first line of 6.25 respectively looks like the 
product of NMHV 1-loop quad-cut x MHV 1-loop quad-cut and MHV 1-loop quad- 
cuts x MHV 1-loop quad-cuts. Therefore, one might think that the right integrals 
to start with look like the product of pentagons that appear in MHV and NMHV 1- 
loop amplitudes. This strategy worked perfectly in the MHV 2-loop case, where the 
amplitude was literary made from double-pentagons whose origin was in the product 
of two MHV-like pentagons. So the natural objects to consider here are the same 
double-pentagons as in MHV 2-loop case and also other double-pentagons that look 
like NMHV 1-loop x MHV 1-loop: 



The numerators of the first three graphs have the same structure as the ones 
that appear in the NMHV 1-loop integrand. We provide the complete expressions in 
appendix B. 

Note that first three diagrams are really represented just by single diagram with 
permuted indices. For instance, the second one can be obtained from the first one 
if we require k > i. So, it is non-planar version of the first graph in the same sense 
as we saw in the last subsection in the case of the log of MHV amplitude. We see 
that these four graphs are in one-to-one correspondence with the first four colored 
graphs in 6.25. If we cut all propagators except (ABCD) we get not only the same 
cuts as are in these colored graphs, but also the support on the correct Schubert 
problems. These integrals are definitely the right ones to start with. In order to get 
the correct field theory answer we have to multiply them by the leading singularities 
of corresponding octa-cuts which are 

Now summing over all allowed indices we get, 
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i i j k 




(6.25) 



where the first diagram really represents three as we mentioned earlier, namely, the 
complete set of cyclically ordered figures 



i I j I k I 




m m i m 



The rest of the story proceeds in the by now familiar way. Simply carrying out 
the sum over the range of indices corresponding to the colored graphs does not give 
the right answer, however, a judicious choice for the range of summation adds the 
correct "boundary terms" to give exactly the right answer, and we finally obtain: 



j I j k 




These two terms represent the general 2-loop NMHV amplitude for any number 
of external particles. The explicit forms of the integrals in term of momentum- 
twistors are presented in appendix B. 
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6.5 3-loop MHV amplitudes 



Finally, we present the integrand for all 3-loop MHV amplitudes. These amplitudes 
were studied in the past, the 4pt formula for the integrand was given in [44] and the 
5pt in [45]. The 4pt and 5pt amplitudes were also determined using BCFW recursion 
and translated into pure momentm-twistor integrals in [1]. However once again our 
new strategy will both yield vastly simpler expressions for these integrands and also 
generalize to all n. 

We begin as always by drawing the colored graphs that contribute to general 
3-loop amplitude. While there are a large number of them, our experience with the 
2-loop NMHV calculation tells us that for the purpose of "translating" the graphs into 
the integrals, one needs to focus on the colored graphs without internal propagators. 
There are just two of these: 




(6.27) 



The colored graphs suggest that the correct 3-loop integral must correspond to "glu- 
ing" together three 1-loop MHV integrals. But these can not be just pentagons 
because of number of internal propagators, we would also need hexagons. Fortu- 
nately, in the "polytope picture" of [9], the most natural form of MHV amplitude is 
written using hexagons. We leave the detailed exploration of this gluing procedure 
to future work. It suffices to say that we can indeed find objects which have support 
on the correct leading dodecacuts (6.27). Having identified them, the magic happens 
again: to get the full 3-loop amplitude, we need only to identify the correct ranges 
for the summations involved. As a result, we can write the general 3-loop MHV 
amplitude for any number of external particles as a sum of two structures, 
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The explicit formulas for these graphs with all numerator factors are given in the 
appendix C. 
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A. Residue Computations in Momentum-Twistor Space 



In section 3, we gave a heuristic argument for the form of the Jacobian in the com- 
putation of the residue of a pentagon integral. The actual computation is essentially 
trivial but it might serve as yet one more way to get used to momentum twistors. 
This is why we carry it out in detail in this appendix. 

Recall that the non-vanishing residue of the pentagon integral for a contour 
which 'encircles' the isolated pole (AB) = (24) is computed using 

r (AB 13) (12 45) (23 45) 

/ (AB 12} (AB 23} (AB 34} (AB 45} (AB 51}' 1 ' ' 

|(AB)-(24)|=6 

As with all multidimensional residues, the entire computation amounts to Jaco- 
bians. Let us choose to expand Za and Z B using the twistors {Z 5 , Z ± , Z 2 , Z4} as a 
basis: this parameterization introduces a Jacobian J = (512 4) -2 . Exploit- 

A£->(512 4) 

ing the GL(2)-redundancy of the integrand, may therefore parameterize Za and Zb 
according to 

Za =Oi\Z 5 + a 2 Z 1 + Z 2 ; 

(A.2) 

Z B =ftZ 5 + /3 2 Z, + Z 4 ; 

Of course, the contour being evaluated corresponds to the choice of maps /j given by 
/ = {(AB 12), (AB23}, (AB 34), (AB 45)}; using these coordinates for Z A ,Z B , the 
contour will be evaluated around the pole at the origin: on = Pi = 0. 

With this, the integral in question has become fully gauge-fixed and concrete: 

M , (Ai? 13)(12 45)3(23 45) 

1 Pl (AB 12) (AB 23} (AB 34} (AB 45) (ABbl) ' 1 ' 

\ai\=\Pi\=e 

Because the contour encircles the origin, the Jacobian appearing the definition of a 
multidimensional residue will be evaluated at the origin. This means that for our 
purposes, we need only compute the maps / to linear-order in (a,, /%) to compute 
the residue. 

Doing this in complete detail, we see that 

/1 = (AB 12) = ax(5412) +... 
h = (AB23} = «i(5413) +a 2 (1423) +... 
f 3 =(AB34}= ^i(2534) +/3 2 (2 134) +... 

U=(AB 45)= /3 2 (2145) +... 
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where '. . .' stands for terms quadratic in a i: fa. From this, it is trivial to read-off the 
Jacobian: 

J = (541 2)(1 423)(2534)(2f 45} = (f 245) 2 (f 234)(2345); (A.4) 

combining this with the rest of the integrand — e.g. (AB 13) / (AB 51) evaluated on 
(AB) = (24)— we find that 

(AB 13) (12 45) (23 45) (2 4 1 3) (1 2 4 5) 3 (2 34 5) 



(AB 12) (AB 23) (AB 34) (AB 45) (AB 51) 

|(AB)-(24)|=6 



(1245) 3 (1234)(2345) 



1. 



B. Explicit form of the 2-loop NMHV Amplitude 

In this appendix, we will provide all the details that go into the formula for the 
n-point 2-loop NMHV amplitude, which can be graphically represented as follows: 



A' 



2— loop 
NMHV 



i<j<l<m<k<i 
i<j<k<l<m<i 
i<l<m<j <k<i 




+ 




i<j<k<Ki 



X + l,k,k + l] 



l NMHV 

4 tree 



+^NMHv('i ■ ■ ■ ij) 
+-^NMHv(^' ' • ' ' 



(B.l) 

Of these two terms, only the first requires any comment, because the second sum- 
mand involves only the familiar double-pentagons which generate the MHV two-loop 
amplitude's integrand, 




(B.2) 



(AB (i-1 ii+1) (j-1 J j+1)) (CD (k-1 k fc+1) fl (/-l U+l))(i j k I) 



(AB i-1 i) (AB ii+1) (AB j-1 j) (AB j j+l)(AB CD) (CD k-1 k) (CD kk+1) (CD l-l I) (CD ll+l) 

(AB,CD) 
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As indicated by the ranges of the summation, the first sum actually represents a 
sum over three distinct cyclic orderings of the labels (i, j,k,l,m), corresponding to 
each of the following cyclically-ordered integrands, 




For each range of indices, there are boundary-terms for which the general integrand's 
numerator must change slightly; these have been indicated in the table above. Given 
the ranges and boundaries indicated above, the numerators for these contributions 
to the 2-loop NMHV amplitude are given by, 

term numerator 



non-boundary (AB i i+l) f] C^i,j,k)) 

A boundary (AB i+l(i-li)f](^i,j,k)) (B.3) 

B boundary (AB i-1 (i i+l) f] (E iJ)fe )> 

A&B boundary (ABi+li-l)(i£j J)fe ) 

where in all these cases the special plane j >k is given by the same object encountered 
at one-loop, but with the arbitrary bitwistor X replaced by (Im), 



_ 1 



(jj+i) ((i k k+i) n (h) - ( k k + l ) ((ij n (h)] • (b.4) 
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C. Explicit form of the 3-Loop MHV Amplitude 



In this appendix, we present the explicit form of the n-point 3-loop MHV amplitude, 
which we represent graphically graphically represented as follows: 



a 3— loop 
•^MHV 




H<«2<il< 
<72<fcl<fe2<«l 



As described in the body of the paper, the 'boundary terms' of the summands 
above require some comment. We will discuss the two topologies separately, starting 
with with the first summand in the equation above. Because when any two of the 
indices become identified in the first graph the wavy-line numerators become ill- 
defined, special consideration must be made for each of the degenerations allowed in 
the range of the summand — that is, all the cases where two or more of the indices 
are identified. Separating each type of such degeneration that is allowed in the first 
summand, 




1 x 



i 



3 Z-^i 

i\<ii<ji< 
<32<ki<k2<h 



Zl(il,i2,jl,j2,h,k 2 ) 



all indices \ 
distinct / 



3x§ E 



l£(ii,i2,ji,j2,k) (k 1 = k 2 = k) 



n<«2<ii< 



3 x| 



1 x 



^l<^2<3^< 
<32<k<h 

E 

il<i2<j<k<ii 

1 E 

i<j<k<i 



l£{ii,i2,j, k) 



i,3, k) 




Here, the overall factor of £ |' reflects the Z 3 -symmetry of the loop integrand (recall 
that every term in the sum is understood to be fully- symmetrized with respect to the 
3! permutations of the loop-variable labels)] although every term in the summand 
has the same factor of |, the boundary terms for which e.g. k\ = k 2 in the sum 
are equivalent to those where ji = j 2 or i\ = i 2 , allowing us to represent all three 



degenerations with a single integrand — T% in this case, and similarly for X:. 



A 
3 ' 
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Let us now carefully define the contributions to this class of graph each in turn. 
First, we have the generic integrand: 



Zi(ii,i2,ji,j2,h,k 2 

for ii<i'2<ji<j'2<ki<k2<ii 




Numerator 
Tr [(n \AB\ i 2 )(ji \CD\j 2 )(h \EF\ k 2 )\ 



Here, we have left implicit the twelve propagators shown in the figure by solid 
lines, and the three 'wavy-line' numerators (AB (ii-liiii+l)f)(i 2 -li 2 i 2 +l)) etc. 
Observe that we have introduced a new notation for remaining tensor components of 
the numerator for this integrand. Letting '•' denote an arbitrary bitwistor, we may 
define a 'trace' over a pair of such auxiliary bitwistors: TV [(a6«)(«cd)] = {abed); 
that is, the trace is nothing but the completely-antisymmetric contraction of bitwistors 
which are dual to a pair of auxiliary bitwistors, which are indicated by '•' in the cor- 
responding formula. 14 

It may be helpful to illustrate the meaning of this numerator using the familiar 
notation of Wick contraction; in this notation, the tensor numerator of X^(ii, i 2 , j'1,,7'2, ki, k 2 ) 
corresponds to: „ 



Tr [(h \AB\ i 2 )( 3l \CD\j 2 )(k 1 \EF\ k 2 )\ = (AB (i x •) n (i 2 •))(CD (ji •) n (j 2 *))(EF (h •) fl (k 2 » )); 

alternatively, the numerator can be written in any one of the following equivalent 
forms (the equality of which offering further justification for calling this a 'trace'): 

Ti[{i l \AB\i2){ji\CD\j 2 ){k 1 \EF\k 2 )} 

= (i2 3i[(j2 h ( (k 2 hA) n f1 (DC)] B)-(A^B); 

= (32 h[(k 2 h ( (i 23l C) fl (SA))) f1 (FE)] D)-(C^D); 
= (k 2 k \U2j1 ( (32 hE) fl (DC))) fl (BA)] F)-(E^F). 



As we will see presently, this numerator will change only very slightly for the 
boundary terms in the summand. Always leaving the propagators and wavy-line 

14 The idea of 'tracing' over auxiliary bitwistors turns out to be a very powerful generalization of 
the four-bracket. Indeed, all the four-brackets in this paper could be translated directly into traces, 
and often with considerable simplification. 
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implicit from the the corresponding figures, the remaining integrands are defined 
according to the following: 



I 2 (ii,i2,ji,j2,k) 

for ii<i2<ji<j2<k<ii 



Z£{ii,i2,3, k) 

for i\<i2<j<k<ii 



• k) 

for i<j<k<i 




Numerator 
TrKh \AB\i 2 )(j 1 \CD\ j 2 )(k \k-l k+l\ k)] 



Numerator 
Tr [(h \AB\ i 2 )(j \j-lj+l\j)(k \k-l k+l\ k)] 



Numerator 

Tr[(i \i-li+l\i)(J \j-lj+l\j)(k\k-lk+l\ k)} 



For the second topology, the boundary terms in the summand lead to just three 
separate contributions that must be specifically addressed. 




h<ji<ki< 
<k 2 <j 2 <i2<h 



all indices \ 
distinct / 



h<ji<ki< 
<k2<j2<l2<il 



2 x l S I 2 B ( i i,ji,k 1 ,k 2 ,i 2 ) (k 2 =j 2 = k 2 ) 



il<jl<ki< 
<k2<t2<il 



lx l Yl l£(ii,ki,k 2 , 



*2 



«i<fci< 

<k2<i2<i\ 



n = 31 = n 
k 2 = j 2 = k 2 



As above, the overall factor of '|' reflects the Z 2 -symmetry of the integrand (we 
remind the reader that each term in the summand is to be fully-symmetrized with 
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respect to the 3! permutations of the loop variables). As before, we have exploited 
the symmetry of the integrand to identify various boundary terms: the degenerations 
h = ji and k 2 = j 2 , being equivalent in the cyclic sum, they can be combined into 
the single summand I 2 — which explains its relative factor of 2. 

With this, we can directly present the three classes of integrands of the second 
topology which contribute to the 3-loop MHV amplitude: 



for U<Jl<fcl<fc2<j2<«2<U 



li(ii,ji,ki,k 2 ,i 2 ) 

for i\<j\<k\<k2<i2<ii 




Numerator 
(AB (i 2 iu' 2 )n0'i-lju'i+l)) 
x(AB (h-lj 2 j2+l)f)(ji fci fc 2 )> 



Numerator 
(AB (i 2 iik 2 )r\(ji-ljiji+l)) 
x(k 2 +lj! h k 2 ) 
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